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Abstract 

We propose and solve a simple but very general quantum model of an SU(2) 
spin interacting with a large external system with N states. The coupling is de- 
scribed by a random hamiltonian in a new general gaussian SU(2)xU(iV) random 
matrix ensemble, that we introduce in this paper. We solve the model in the large N 
limit, for any value of the spin j and for any choice of the coupling matrix element 
distributions in the different possible angular momentum channels I (and provided 
that the internal dynamics of the spin is slow). Besides its mathematical interest 
as a non-trivial random matrix model, it allows to study and illustrate in a simple 
framework various phenomena: the decoherence dynamics, the conditions of emer- 
gence of the classical phase space for the spin, the properties quantum diffusion in 
phase space. The large time evolution for the spin is shown to be non-Markovian in 
general, the Markov property emerging in some specific case for the dynamics and 
the initial conditions. 
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1 Introduction 



In this paper we introduce and solve a simple but general model of decoherence for a 
very simple system: an SU(2) quantum spin coupled in a generic way to a large external 
system, which plays the role of the environment (external bath, other microscopic degrees 
of freedom, etc.). The purpose of this model is to provide a simple but quite general and 
exactly solvable model, in order to discuss and illustrate in a pedagogical way some of the 
basic aspects of decoherence and of the emergence of classical degrees of freedom in a sim- 
ple quantum system. There is of course already a very large scientific literature on these 
problems and this study shall not bring revolutionary insights on these subjects. There are 
numerous excellent reviews or textbooks on decoherence |Schlosshauer, 20 07||Zurek, 2003 



|Joos et al., 2 003 1, and more generally on irreversibility and dissipation in open quantum 
systems | Weiss, 2008||B reuer and Petruccione, 2006] . Many studies have focused on sim- 
ple but realistic physical systems like a single free massive particle (Brownian walk) or 
an harmonic oscillator. A quantum spin is in fact simpler than a particle, since its phase 
space is compact and since the SU(2) symmetry induces many simplifications. The spin 
model presented here does not aim at being a realistic one for some specific physical sys- 
tem in relation with experiments. It incorporates features already introduced and studied 
by previous authors: quantum spin, coherent states representation, random matrix en- 
sembles. But it also presents novel features and allows to obtain new and exact results. 

Firstly our model describes a general spin, and the results that we obtain are valid for 
any value of the spin j, allowing to study the whole range going from the case j = 1/2 
(the most studied case of the two-level system, i.e. of the q-bit) to j = oo (the classical 
limit where the spin becomes a classical top). The study of dissipation and decoherence 
in spin systems goes back to |Shibata and Saito, 1975[ Takahashi and Shibata, 1975| , but 
since then most studies have focused on the spin 1/2 case, i.e. of the two level system, 
both for simplicity and for the obvious connection with experiments. See however Sec. 
VII of Si run/ et al., 2002 for a discussion of decoherence for a large spin j. 



Secondly we try to describe the interaction between the spin S and the external sys- 
tem £ in the most general situation, by considering an interaction Hamiltonian belong- 
ing to a random matrix ensemble. The idea of using random matrix theory (RMT) 
for such problems is of course not new. It can be found already in |Mello et al., 1988, 



Lutz and Weidenmiiller, 1999| . However, as already stated, most studies involve simpler 
systems (for instance the two state system corresponding to j = 1/2), and some simple 
random matrix ensembles like a single GUE ensemble or au Gaussian ensemble of band 
matrices. In addition in almost all studies, the interaction Hamilonian is taken to be of 
the form H mt = Us <8> Vg, Us being a well chosen operator for S (the coupling agent), 
and Vg the random Hamiltonian for £. Here we look for the most general random matrix 
model which can describe the coupling of an SU(2) spin to an external N states sys- 
tem, irrespective of the value of the spin j. This leads us to ask the question: what are 
the most general random matrix ensembles which are invariant under the global unitary 
group U(iV) for the external system £, and the SU(2) symmetry group for the spin S. We 
solve this problem and define a class of Gaussian SU(2)xU(iV) matrix ensemble, that we 
denote GL^xatE. This class of random matrix ensembles is new (as far as we know) and 
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has interesting properties. In particular they can be formulated and visualized in terms 
of random Wigner and Husimi distributions on the Riemann sphere. 

Thirdly, in the limit iV — > oo (i.e. when the external system becomes large) we are 
able to write closed equations, and to write in a fully closed form the evolution func- 
tional for the density matrix ps{t) of the spin S (once the trace over the external system 
£ has been taken). In other word we do not need to make any approximation and to 
write any master equation for the evolution of the density matrix. Our calculations rely 
on the classical methods developed in the physics litterature in |Brezin and Zee, 1994| 
|Bre zin et al., 1 995||Zee, 1996 1 to study sums and products of random matrices, and are 
closely related to the techniques of free probabilities in mathematics (see |Voiculescu et al.~ 992 



and references therein). They are also in fact closely related to the methods used by 
|Me llo et al., 198 8 and more recently by |Lebowitz and Pastur, 2004| |Lebowitz et al., 2007| 
for studying analytically the dissipation in two level systems. 

One important limitation of the model discussed here is that the internal dynamics 
of the quantum spin S is completely neglected. This means that our results apply in 
the specific case where the internal dynamics of S is much slower than the dynamics of 
the external system £ and than the dynamics on the whole system S+£ induced by the 
coupling. We are thus studying the case of "interaction dominated decoherence" already 
put forward by | Str unz et al. , 2002] . This case is relevant for the study of decoherence 
and of quantum diffusion effects, but not for dissipation. 

Rather than discussing the whole literature in this introduction (this would make it 
much too long), we prefer to present our model now. Our results will be discussed and 



compared with previous ones at the end of the paper in section 6.2 (the interested reader 
can go directly to this section). The paper is organized as follows. 

In section [2] we introduce our model, define the random matrix ensembles that we shall 
use, and discuss the connection of these ensembles with the theory of Wigner and Husimi 
representation for spin operators. 

In section [3] we study the evolution functional for the spin density matrix, for a general 
choice of random coupling Hamiltonian, characterized by a set of variances A(/) for the 
couplings in the different angular momentum channels I. We first write recursion relations 
for the evolution functional and the associated resolvent operators in the large N limit 



in section 3.2 We show that these relations takes a simple closed form in each I chan- 



nel when expressed in angular momentum components (their Wigner transform double 



Fourier components). The general solution is given in section 3.3 It involves a universal 
decoherence function M(t,Z(l)) which depends on the time t and of a parameter Z(l) 
which is some linear transform of the variance distribution A(Z'), involving the SU(2) spin 
structure through the Racah 6j-coefficients. The structure of this decoherence function 
as a function of the parameters A(/) and various scaling limits are discussed at length in 



section 3J3 and illustrated in Appendix [Aj The relation between the parameters of the 



decoherence function and the norm of operators built out from the Hamiltonians and the 



spin operators is discussed in section 3.6 



In section [4] we study explicitely the decoherence of spin states in our model. We 
concentrate on the case where the dynamics of the external system £ is also given by a 
GUE Hamiltonian (for simplicity) and when the initial state for £ is generic. For large 
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values of the spin j this is of course naturally discussed in term of spin coherent states, 



whose properties are recalled in section 4.1 In section 4.2 we compute the various time 



scales for the system (which corresponds to Td yn . a dynamical timescale of the whole 
system, Td ec . a decoherence time scale and Tdi SS . a dissipation time scale) as a function of 
the parameters of the model. We then compare the evolution of generic non coherent 
states (random states), coherent states and simple superpositions of coherent states (cat 
states). Our solution shows explicitly under which conditions (choices of the A(Z)) the spin 
coherent states emerge as robust states under decoherence, and therefore are the pointer 
states which allow to describe the spin system and its dynamics in semiclassical terms. 
Our explicit solutions allow very easily to illustrate the dynamics of decoherence through 
the dynamics of the Wigner representation of the spin quantum states. This is done in 
section 4J3 with several pictures and animations. In section 4A we study the large time 
evolution of the semiclassical coherent spin states. We show that through the coupling 
with £, they undergo a slow quantum diffusion process on the classical phase space (the 
sphere <S 2 )- We show that this diffusion process exhibits universal self similar properties, 
but is non Markovian at all time scales. Hence in our case the evolution functional cannot 
be approximated by a master equation local in time. 

In section [5] we extend our calculations to the more general case where the dynamics 
of the external system £ is given by an arbitrary Hamiltonian with a continuum energy 
spectrum, and when the initial state for £ is not taken to be a random state, but is 
an energy eigenstate \E). The coupling between the spin and £ is still described by a 
GU2xAfE Hamiltonian. In section 5.2 we show that we still can compute explicitly the 
influence functional, but that it depends on more parameters, in particular on the energy 
E on the initial external state. The discussion of decoherence in the simple case when 
the spectrum density p(E) for 8 is a semi circle (the explicit results are simpler in this 
case) is done in section 5.3 Finally in section 5.4 we concentrate on the case where the 
internal dynamics of external system £ is much faster than the dynamics induced by the 
coupling between the spin S and the 8. In this case we show explicitly that the large 
time evolution of coherent states is still described by a quantum diffusion process, with 
depends on the initial state for £ we start from. The quantum diffusion process has the 
property that if we start from an energy eigenstate \E) for £ (times a coherent state for S 
of course), the quantum diffusion process is now Markovian, and coincides with a diffusion 
process local in time and space (i.e. the Wiener process). The effective diffusion constant 
in phase space depends on E, and takes a simple form related to the Fermi golden rule, 
indicating a fluctuation dissipation relation. When we start from a general state for 
£ , the diffusion process is not Markovian, but is argued to be described by a randomized 
Markov process. We propose a simple explanation for this randomization, through the 
decoherence for the energy states of the environment induced by the now large spin. 

In section [6] we first summarize the main results on our model, then put it in the 
context of the existing literature (section 6.2), and finally discuss open problems and 
possible generalizations of our model (section 6.3). 
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2 The model 



2.1 A spin coupled to its environment 

We start from the standard paradigmatic system: a small system S coupled to a large 
system £ (the environment). We take for S a single SU(2) spin with value j. This 
has several advantages, we can use the coherent states representation and the techniques 
of group representation theory, and we can easlily study the emergence of the classical 
degrees of freedom in the large spin limit j — > oo. The coupling between S and £ is 
formulated and studied with random matrix models. The advantages are that we can 
use the mathematical tools of random matrix theory and that the interactions that we 
consider have are generic. 

The Hilbert spaces for the subsystems S and £, and the whole system S + £ are of 
course 

H s = C 2]+l , H S = C N , U = U s ®U e (1) 

We take as an orthonormal state basis for Hs the S 3 eigenstates {\m), m = —j, ■ ■ ■ , j}, 
the S^jfi = 1, 2, 3 are the spin operators (the generators of su(2)). An orthonormal basis 
for Hs is taken to be {\a), a = 1,N}. The Hamiltonian of the whole system is written 
(altough this separation is partly arbitrary) as a sum of three terms 

H = H s ® Is + H ss + 1 5 ® H e (2) 

Hs and describe the internal dynamics of S and £ respectively. H$s describes the 
coupling between S and £ and is usually taken to be small (fast internal dynamics). 

We are mostly interested in the decoherence of spin and we therefore assume the the 
dynamics of the spin in slow compared to the other dynamics. This important simplifi- 
cation will be discussed later. Therefore we simply set 

H s = (3) 

Moreover we first assume that the internal dynamics of £ does not play a special role and 
is at most as fast as the dynamics of the coupling S + £. We shall therefore recast Hs 
into Hss, thus setting simply 

H g = (4) 

We shall go beyond this simplification in section [5] In the rest of this section and in the 
following section [4] we shall simply denote Hss by H int or simply by H. 

2.2 Gaussian SU(2) xU(iV) random matrix ensembles 

We are interested in the generic situation where H int is taken to be a typical (or generic) 
Hamiltonian allowing the emergence of the spin as a classical degree of freedom. Therefore 
we assume that H mt is a random Hamiltonian belonging to some M x M random matrix 
ensemble, with 

M = (2j + l)N . (5) 
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As explained in the introduction, in most models H\ nt is taken to be of the form 

#int = U s ®V e (6) 

Us being an appropriately chosen operator for S (the coupling agent), and Ve a random 
Hamiltonian for S (sometimes it is a sum of a few terms like this). Here we look for 
a general random Hamiltonian, simply assuming (mostly for simplicity, but there are 
some physical motivations too in some specific models)) that it belongs to some Gaussian 
ensemble. The question is thus: can we caracterise the most general Gaussian ensembles 
(random matrix probability distributions) which are invariant under: (1) the SU(2) group 
acting on %Si (2) generic unitary transformations U(iV) acting on Tie*? This does not 
mean that H is invariant, simply that the distribution of the H's is invariant. 



2.2.1 Spin decomposition of operators 

To answer this question, we simply use standard tools of group representation theory. 
We refer to |Varilly and Gracia-Bondia, 1989| for a clear reference to the theory of SU(2) 
representations, coherent states and the theory of Wigner functions, that we shall use 
heavily . Let us first concentrate on the spin part, i.e. on the subspace %$■ Any operator 
A acting on Hs, i-e. A G B(%s) — A^2j+i((C) can be considered as an element of the 
tensor product (or Kronecker product) of two spin j representations of SU(2). Of course 
such a product is a sum of the irreducible representations 



j®j =0©1 



2j 



2j 

0> 



1=0 



Hence the operator A can be decomposed into its spin I components A"> 



(7) 



■2.) 



(8) 



The matrix elements of the A^ are 



42 = (r\A m \ s ) = J2(-iy 



21 + 1 1 2 I 



m=-j 



2j + 1 \s —m 



(9) 



where the coefficients W^' m ^ are given in terms of the matrix elements A rs = (r\A\s) of 
A by the inverse transform 



W 



(Lm) 



E 



r,s=-j 



21 + 1 /j I 



2j + 1 \r m 



J 



(10) 



The 



j I 
mi m 2 



m 3 



are the SU(2) Clebsch-Gordan coefficients. Of course only the single 



terms such that m = s — r contribute in the sums. 
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The spin / component of A satisfies 

3 

[S, [S, A^}} = J2i S ^ [ s "> = l ( l + ( U ) 
The coefficients are complex but satisfy the conjugation constraint 

w V,m) = ( _ 1)m ^- m ) ( 12 ) 

2.2.2 Relation with the Wigner and the Husimi distributions 

Wigner representation: The coefficients W A l ' m ^ are nothing but the (l,m) coefficients 
of the Wigner distribution W^n) associated to the operator A in the basis of spherical 
harmonics Y l m (n) on the unit sphere <S 2 . The Wigner distribution is 

Wa(u) < m) Yi m (n) (13) 

1=0 m=—j 

With the normalisation for the Y™ 

[ d 2 nYr{n)Yf 7 {n)=5 hll 5 m ^ , d 2 n = d9 d<P sin(0) (14) 
Js 2 

With these normalisations we have 

ti(AB ] ) = j d 2 n W A (n) W^ftj = W%' m) W%' m) (15) 

^ 2 l,m 

and for the unit operator and the trace 

Wf m > = V2j + l Mm,o , tr(A) = v^JTT <' 0) (16) 
NB: In the litterature the normalization for the W A is often such that 

tr(AB^) = 2 ^-^J d 2 n W A (n) W B (n) (17) 
It corresponds to the change of normalisation in the definition of the Wigner distribution 



W A (n) = s /^ Y W A (n) (18) 
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Husimi representation: Note that we can also construct the Q-symbol, or Husimi 
representation 

<w*o = £ wf m) *r(tf) (19) 



l,m 

and the P-symbol 



;•) W^' m) K"i;n (20) 



They are such that 

tr(AB t )= / d 2 n Q A (n)P^) (21) 

The Husimi function Qa(^) corresponds to the "physical" probability distribution of spin 
in phase space, since with the proper normalisation 



/ 4tt 

(n\A\n) = Qa(^) , Qa(^) = y 2 • _|_ ]_ Qa(^) , |^) coherent state (22) 
2.2.3 Gaussian SU 2 random matrix ensembles: 

In order to construct the most general Gaussian SU(2) invariant ensemble on the self 
adjoint matrices in M 2 j+i(C), we simply have to take separately each matrix component 
A® as Gaussian independent random (non-commutative) variables with zero mean and 
variance depending on I. Namely, we take the W^ 1 ' 171 ^ to be independent Gaussian variables, 
subjected only to the Hermiticity constraint 

w (i,m) = -"») (23) 

which ensures that the A are Hermitian operators, and to the fact that the variance 
depends on / but not on m, which ensures SU(2) invariance of the distribution. More 
precisely, we take the Ws to be for m = 

W m = A {1) (24) 

and for m > 

W (l,m) = B (l,m) + iC (l,m) ^ ^(I-m) = (_-Qm ^(I,m) _ iC (l,m)j for Q < m < / (25) 

and to take for the A®, B^ l ' m ^ and C* 7 '" 1 ) random Gaussian independent variables with 
zero mean and mean square extend A(Z) depending only on Z. The cumulants are 

EL4 W ] = E[B {l > m) } = E[C {l ' rn) ] = (26) 

E[A {l) A {l) ] = E[B (l ' m) B {l ' m) ] = E[C {l ' m) C {l ' m) ] = A(Z) (27) 

All the others cumulants being zero. The A(/)'s are a collection A of 2j + 1 positive 
numbers 

A(Z)>0, Z = 0, - - - 2j (28) 
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which completely characterize the SU(2) gaussian ensemble. This distribution is given by 
the Gaussian probability measure on self adjoint (2j + 1) x (2j + 1) matrices 



£> A [A] oc dA exp 



2j 

E 



2A(Z) 



tr 



(29) 



where dA is the standard flat measure, so that E[F[A|] = J T>& [A] F[A]. 
In this GU2E ensemble, characterised by A, the "propagator" V rs ^ u is 



2.7 



V rsM = E [A rs A tu ] = S s - r!t - u 2J A(/) 



21 + 1 /j I 



1=0 



2 j + 1 \s r - s 



r \t u — t 



(30) 



Of course, if all the A(Z) are equals to the same A, one recovers the standard GUE 
ensemble for (2j + 1) x (2j + 1) matrices, with 



IE [A rs A tu ] G ^ E — A 5 r>u 5 s j 



2.2.4 Gaussian SU^xUat random matrix ensembles: 



(31) 



We now take into account the external system £ and want to characterize the most general 
Gaussian ensemble of self-adjoint matrices H G B(H) = M^+iwCC) which is invariant 
under SU(2) (acting on l-i s as the spin j representation) and under U(iV) (acting on 
as the fundamental representation). 

The solution is simple. We make the same decomposition w.r.t. the spin sector Hs , 
keeping the sector Hs untouched. H is decomposed into 



H 



(i) 



(32) 



The matrix elements of H are now denoted as 

F r f = (ra\H\s(3) , 

and can be written as 



\ra) = \r) <8> |a) 




W. 



(l,m) 



a 3 



The hermiticity constraint H = for H reads now 



W. 



Q,m) 



'-l) m w 



(l,-m) 



(33) 



(34) 



(35) 



The GU2 ) at-E' ensemble is obtained by taking the W^T to be Gaussian independently 



distributed random variables, subjected to the constraint of eq. (35) , with zero mean 



and a variance depending only on /, but not on m, a and (3. This ensures the invariance 
under the group SU(2)x U(iV) of the probability distribution on H. The "propagator" T> 
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is now the product of the propagator T> for the GU2E ensemble (given by eq. (30)) times 
the propagator for the standard GUE model. 



rs.tu 



E 



Aa/3 ajS 
^rs *Hu 



V 



rs,tu $a,S ^/3,7 



(36) 



This ensemble is characterized by the same ensemble A = {A(Z), Z = 0, • • • 2j} of positive 
parameters. 

Let us note that the sector I = gives a Hamiltonian independent of the spin, 
since its matrix elements are of the form 



Hence can be written as 



u r,s rL (g) 



(0) 



(37) 



(38) 



where if( ) is a random Hamiltonian for S, whose distribution is given by the GUE 
ensemble with (variance) 2 = A(0). 

Finally we note that the Hamiltonian can of course be rewritten in the general form 



H 



(l,m) 



(l,m) 

s 



(l,m) 



(39) 



where the Dg' m ' are some fixed spin operators (related to the Wigner D-matrices) and 
the Wf m) some random operators on S. Previous studies of spin decoherence have dealt 
with simpler interaction Hamiltonians Use with only one D ® W term or a few (typically 
3 when dealing with a / = 1 interaction). Here we keep all the possibles terms in the 
decomposition 39 But we shall stay with the explicit form 34 for the decomposition of 
H. 



2.3 Rescaled distributions with TV and j 

We thus take for our model a single spin j, with for simplicity no internal dynamics H$ = 
and for Hamiltonian H$e describing the coupling between the spin and its environment 
a random Hamiltonian in the SU(2)xU(A r ) random Gaussian ensemble characterized by 
the family of coupling amplitudes A = {A(Z), / = 0, • • • 2j} in the different spin channels 
/. In this paper we are interested in the limit iV — > 00 (large environment) and j — > 00 
(classical spin). In these limits it is adequate to rescale these amplitudes (i.e. the time 
scale for the evolution of the system). For clarity we define here these rescaled parameters 
as 

A(/)=iVA(/) , A(/) = ^ T A(/) (40) 

As we shall see later, the large A^ limit is obtained when keeping the A(Z) of order 0(1), 
while the large A^ and large j limit is obtained when keeping the A(/) of order 0(1) 
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3 The evolution functional 



3.1 General framework 

We start at time t = from a separable quantum state 

p(0) = /(0)®/(0) (41) 

We do not specify at that stage if the system S is in a pure or in a mixed quantum state. 
Since the ensemble for the Hamiltonian H is chosen to be invariant under arbitrary unitary 
transformations U G U(iV) acting on TLg, and since we shall compute only averages w.r.t. 
the distributions of H, the precise initial quantum state p £ (0) is not important for what 
we are interested in. We may choose it to be the maximal entropy state 

/(0) = jjU (42) 

After evolution from time to t > of the whole system under the time independent 
Hamiltonian H, the reduced density matrix for the system S is 

p s (t) = ltr £ (e- itH (p s (0) ® l £ y m ) (43) 

We are interested in the evolution of the model for a "typical" Hamiltonian H in our 
GU2x7v-E ensemble, and we shall therefore make averages over H in this ensemble. The 
average reduced density matrix 

^(t)=E[p s (t)] (44) 

is sufficient to study the observables of the subsystem. However in order to study deco- 
herence and the evolution of the intrication between the systems S and £ , we must study 
functions of the reduced density matrix such as the von Neumann entropy S or the purity 
P (or whichever quantities you prefer) 

S = -tx s [ P S \ogp s } , P = ti s [p s p s ] (45) 

It is well known that simplifications occur for such observables in the limit iV — > oo 
where the number of degrees of freedom of the external system £ becomes infinite, with 
a adequate rescaling of time and of the couplings. Indeed, as we shall discuss in the next 
section, for "reasonable" functions F of the density matrix p s (and at least for polynomial 
functions), in the large N limit we have the factorisation of the expectation of products of 
traces with respect to £ due to the fact that in this limit, only planar diagrams contribute 
in the perturbative expansions. This means that we have in fact Q 

FG^) =E[F(p s )] = (1 + 0(N- 2 )) (46) 

1 In the rest of this paper we shall rather use the physics notation F for the average of F rather that 
the probabilistic notation K[F] for the expectation of F, but this is of course the same quantity. 
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This phenomenon is known in physics as "factorisation" or emergence of a large N master 
field. In mathematics it is known as the phenomenon of "concentration of measures". It 
leads to the formulation of the N — » oo limit in terms of free probabilities. 

It turns out that these kind of Gaussian ensembles are self-averaging (this is discussed 
at length in the mathematical literature, see for instance |Lebowitz et al., 2007|), but we 
shall not elaborate further this point. 

Thus we have to consider the evolution operator Ai (t) for the density matrix 

7ru(t) = J2 M ruAt)pU0) (47) 

st 

Ai(t) corresponds to the Feynman- Vernon influence functional, it is of course a strictly 
positive linear trace preserving application on M 2 j +1 (C), thus it is also called a quantum 
channel, or a POVM. Its matrix elements are given by 

MruAt) = ^^(ra|e-^| S /3)(t/3|e^| M a) (48) 

3.2 Perturbation theory and the planar limit 
3.2.1 Resolvents 

It is more convenient to consider the double resolvent 

g rUiSt (x, y) = ^2(ra\{x - Hy l \sf3){tf3\(y - Hy l \ta) (49) 

Formally (integration paths to be discussed later) 

M ru , st (t) =f^ K f^ K e~ it{x - y) QmA^ V) (50) 
We consider first the single resolvent 

% rs{x) = ^^(raKx-i/rV) (51) 

a 

To compute these generating functions, we use the standard diagrammatic techniques 
(see for instance IMello et al., 1 988 1| Zee, 1996 1 and references therein). We represent the 

j_ t-»ch9,70 i ii £ j- v mi liii- j j i j i i. 



propagator X>"f^ as a double fat line. The dashed lines represent the two external tensors 

a 5 



rs,tu s — { 

p y 

Figure 1: Diagramatic representation of the propagator T> 



a/3,75 

rs,tu 
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5 a> s and 8p tl , and indicates that the £ indices (greek letters) are conserved. The black 
ribbon represents the spin tensor T> rSttu - It indicates that the spin indices (roman letters) 
are mixed. But the difference between the left and right indices is conserved. 

s — t = r — u (52) 

3.2.2 Recursion equation for the single resolvent 

To compute T-L TS (x) we expand in a power series in x~ l 

oo 1 

% rs {x) = ^aT 1_fc — J2( ra \ Hk \ sa ) ( 53 ) 

k=0 a 

and use Wick theorem to compute the average H k . We get a sum of contributions associ- 
ated to diagrams of the form depicted in Figj2j The propagators form arches above a line 
going from r to s. Each arch gives a term T> proportional to the A(/), and each closed 




Figure 2: The first diagrams for % rs (x) 



dashed loop gives a factor iV (sum over the £ indices a). So each diagram is of order 



x 



-l-2#arches ^j#arches jy#loops-l _ x ~l ^ x -2j^j jyj #arches ^ 



(54) 



X being the Euler characteristic of the fat diagram. [A] means any A(/). Thus in the large 
iV limit, only planar diagrams (rainbow like) survive, provided we rescale the variances 
by N 

A(0 = JVA(0 = 0(1) N^oo 
% rs {x) then satisfies the recursion equation 



U 



(55) 
(56) 



t.u.v 



with 



2.7 



V 



rt,uv 



NV rt ,uv = 8s-r,t-u / ] A(l) 



21 + 1 /j I 



1=0 



2j + 1 \s r - s 



Ah i 

r/ \J u-t 



(57) 



This recursion equation is depicted graphically in Fig. [3} The solution is of the form, in 
fact required by SU(2) invariance 



H r ,(x) = 8 r ,H(x) 



(58) 
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Inserting this ansatz 58 into 56 and 57 we obtain the simple recursion equation for TL(x) 



2j 



T-L(x) = x 



-l 



9/ 1 

l U(x) 2 A with A = ^2=f^A(l) 



1=0 



2j + l 



(59) 



Hence 



U{x) = ^ ( x - \/x 2 - 4A ) = x" 1 Cat(A.T" 2 ^ 



where Cat (2) = ^ z n C n is the generating function of the Catalan numbers C n . 



(60) 





Figure 3: Graphical formulation of the recursion relation eq. (56) for T-L rs (x 



3.2.3 Recursion equation for the double resolvent 

Now we can compute the function G r u,st{x,y). It is given by the sum of the planar 
diagrams of the form given in Fig. |4j It thus obeys the recursion equation 




r st u 

Figure 4: The planar diagrams of Q r u,st{x,y) 



Gru.sl ('I': if) ^ ftr.sUlu (?/) + 2* T^rv,wx^-vw ( i %)Gxu,st (•£ j V) + V ^rv,wxGvw,st^ j Z/)^xu(?/) 

(61) 

To solve this equation, it is better to use its SU(2) invariance properties, and to rewrite 
it for its "double Wigner transform" coefficients 



W^> mi ^ h > m *\x, 



y)= E E 



'2Zi + l /j I 



r,u=—j s,t=—j 



2j + 1 \umi 



'2*2 + 1 /j k 

2j + 1 \s m 2 



Gru,st(%, y) 

(62) 
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n 



u r s t 



u r v w x s t 



I I 

r v st w x 




Figure 5: Graphical formulation of the recursion relation eq. (61) for Gru,st(x,y) 



Indeed, we reexpress the initial propagator, V> r s.tu (given by 57) in the (s,t) — > (u,r) 
channel as 



T-^rs,tu T-^ru,st 1.6. 



x: 



y 
A 



(63) 



The corresponding double Wig ner transform of D ru s t is 



W, 



(il,mi),(/ 2 ,m 2 ) 



V 



E E 

T,U=— j S,t = —j 



'2h + l /j I 



2j + 1 \u mi 



j\ 2l 2 + l/j I 



r j W 2j + 1 \s m 2 



^ I T-^ru,st 



(64) 



Using the original expression 30 for T> rstu , this sum is rewritten as a multiple sum over 



products of four Clebsch-Gordan coefficients. After some SU(2) algebra, it can be reduced 
to the simple form 



W 



(h,mi),(h,m2) 



V 



dh,h ^mi+m.2,0 (— l) mi A(/i) 



(65) 



with A(/i) given by 



2j 



A(h) = A(0(2i' + i)(-i) 2i+l ' +h y. j I 



l'=0 



(66) 



eq. (59) 



J4 J5 Je 



where <{ ^ ^ 2 ^ is the Racah 6-j symbol. In particular, A(0) is nothing but the A of 



2.; 



a(o) = A = ^|±1a(/) 



1=0 
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2j + l 



(67) 



Remember that the A's are just the original A's rescaled by a factor of N, A(Z) = NA(l). 

The constraints l\ = I2 and ni\ + m-i = are very important! They just express the 
SU(2) invariance of T> rSttu , i.e. the initial contraint s — t = r — u. But when we take 



the Wigner transform of the recursion equation 61 for Gru,st(x } y), one see that they are 
preserved by the equation and that the solution for Wg 1 ' mi ^^ 2 ' rn2 ^ x , y) must be of the 



same form as 



W^' M (x,y) = 5 h>l2 5 mi+m2>0 (-l) mi & h Xx,y) (68) 



We can now come back to the recursion equation 61_ Using 62 and 68 it factorizes 
into independent equations for each Q^ ll \x,y) 

g {l) (x,y) =x^(i/)+x- 1 A(0)^(x)^ ) (x,i/)+r 1 A(0^ (0 (^y)^(z/) (69) 



The solution is simply (using the explicit form 60 for H) 



gU(x,y)= n{x)n{y) (70) 
1 - A(l)H(x)H(y) 

3.3 Solution for the evolution functional 
3.3.1 General form 

We can now obtain the influence functional M.(t). As we shall see, the functions Q®(x, y) 



are analytic in x and y around 00, and have a cut in the x and y planes along [— 2 v A, — 2 v A] . 
We can integrate in x and y along a closed anticlockwise curve around the cut to obtain 
the double Wigner transform of the influence functional Ai ru ,st(t) 

w (h,m l)M{t) = Kh 6m ^ Q ( _ 1)mi fa {h){t) (n) 

where 

M {l \t) = f ^ f ^e-^-* g«\x,y) (72) 

Therefore, the evolution of the reduced density matrix p s (t) becomes a separate simple 
linear evolution in each (l,m) sector when one considers the components of its Wigner 



transform. More precisely, if W^ m \t) is the (/, m) harmonic, given by 10, we have simply 

W^ m) (t) = M {1) (t)W^ m) (0) (73) 
Then using M we can reconstruct p s (t) in the \r)(s\ basis. 



3.3.2 General decoherence function 

Thanks to the SU(2) invariance, the evolution functional reduces to a single function in 
each I sector. This function depends on time and on the distributions of the A(Z) which 
measure the strength of the coupling between the spin and the environment in the different 
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angular momentum sector I. This function depends in fact only on two parameters, since 
it can be rewritten as 

M®(t) = M(t/T ,Z(l)). (74) 

where r is a time scale 

21 21 + 1 ~ 



r = l/VA(0) , A(0) = £ 27XT A (0 ( 75 ) 

i=o ■> 



and Z(l) a parameter depending on the angular momentum I 



2/ 



A(0 ,to [3 3 I 



J 3 V 



Z{1) = ^ = ^ (76) 

A(0) I mK{r) 

l'=0 

and the decoherence function M(t, Z) is simply 

if (x) is nothing but the resolvent of the standard normalized GUE ensemble. It is the 
Hilbert-Stieltjes transform of the Wigner-Dyson semi-circle density distribution. It has 
a cut along the interval [—2,2], behaves as x~ x at oo and its modulus is < 1 for 

x G C\[— 2, 2]. Hence the function M(t, Z) is well defined for any real t, and analytic in 
the disc \Z\ < 1. We shall discuss its properties below. 

Thus we have a completely closed and simple formula for the evolution functional of 
a spin coupled to a large environment via a random coupling Hamiltonian which belongs 
to an SU(2)xSU(N) invariant ensemble. Our formula is valid for any value of the spin 
j, going from j = 1/2 (the q-bit or two level system) to j — > oo (the classical spin), and 
for any distribution A(/) of the strength of the couplings as a function of the total spin I 
exchanged via the interaction. 

Our result separates in two parts: (1) the universal decoherence function M(t, z) which 
comes from the RMT part of the calculation; (2) the parameters A(Z) which depend 
linearly from the initial distribution A(/) of the couplings as a function of the angular 
momentum I, which come from the SU(2) group theory part of the calculation. They give 
simply the time scale To and the parameter Z(l) 

3.4 Properties of the decoherence function M(t,Z): 
3.4.1 Analytic representation 

Making the standard inversion of variables x — > H as in |Ze e, 1996] (i.e. going from the 
Green function H{x) to the so-called "Blue function" B(w)) 



w = H(x) <^=^> x = B(w) = w + w 1 (7£ 
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we rewrite M as 



M(t, z) = I ^ / e-i^)^) ^!-^- 1 )^-^ 1 ) (7g) 
/ 2i7r J 2m I — ZW1W2 

where integrating along the cut [—2, 2] in W7\ amounts to integrate along the unit circle in 



79. We can use it to obtain the double (t, x) series expansion of M(t,z) which is found 



(after a bit of algebra) 

00 m 

m!(m + l)!(m — n)!(m + n + 2)! 



M(t, z) = V fW(-l)^ „ 2(2m + 1)(n +,^ (2m)! „ (80) 



m=0 n=0 

Thus M(t, z) is a generalized hypergeometric function of the two variables t 2 and z. It is 
depicted on fig. ^ 




Figure 6: The decoherence function M(t, z) as a function of the rescaled time t and the 
parameter z G [—1, 1]. Decoherence is fast when z < 1 but very slow when z ~ 1. 



3.4.2 Small t limit 

The small £ behavior of M(t, z) is 

M(t,z) = l + t 2 {z- 1) + C(t 4 ) (81) 
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3.4.3 Large t limit 

The large time behaviour of M(t, z) is most easily calculated from the integral represen- 



tation 79 by using the steepest descent method at the saddle points wi = ±1, w 2 — ±1- 
We obtain an algebraic decay as t~ 3 , with an oscillatory term negligible when z — > 1 and 
dominant when z — > — 1. 

3.4.4 The z ->• 1 and i(l - z) = 0(1) scaling 
When z = 1, we have in fact 

lim M(t, z) = 1 (83) 

but the function M(t, z) takes a scaling form when z — > 1 while £ is large. In fact 

M(t,z) = ^(t') with *' = *(!-*) in the limit t' = 0(1) , z ->■ 1_ (84) 



Indeed, in this limit, the measure in eq. (79) concentrates around w 2 = W\ and the integral 



becomes a single integral representation. We get a simple hypergeometric function 

dw 1 „, x „ 2 r /2 



/ 2i7T 2 W 1 ^ TT y_ ff/2 



1 /-2 

--- / di^e" 1 '^ (85) 

27T 



-2 

whose series expansion is explicitely 



This function is depicted on fig. ([7]). Its asymptotic behavior is 

1 fi 1 
#(£') = 1 - -JL*' + £>(t' 2 ) £' , V(0 = " t'~ 3 + C(^" 4 ) *' ->• 00 (87) 

3 7T 7T 



3.4.5 The z -»■ 1 and i = 0(1) scaling 

Note that \l/(t) is linear in t at small time, not quadratic in t like M(t,z) for z < 1. For 
z = 1 — e close to 1 (e <C 1) but t of order 1, the function M(t, z) behaves as 

M(t,z) = l + e$(t) + 0(e 2 ) (88) 

with $(£) = 1 — iF 2 (— h 1,2; — 4t 2 ) a universal non-linear function which behaves as 

$(t) = -t 2 + C(t 4 ) when t , $(t) = --i + 0(1) when i -> 00 (89) 

7T 

Therefore, the crossover between the non-linear regime [88] for small t and the linear 
regime 84 for large t ~ e _1 occurs in a domain of t of size 0(1), hence in a very small 



interval in t' of size 0(e) 
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Figure 7: The scaling function \I/(t) for the decoherence function when z — > L 




-2.0 L 



Figure 8: The function $(t) of eq. (88) and eq. (89) that describes the small t and z ~ 1 



(hence I <C j) behavior of the function M(t, z). 



3.5 Z(l) dependence on the choice of couplings A(7) and on the 
total spin j. 

3.5.1 Various examples 

Now we study how the time scale To and the parameters Z(l) (which govern the dynamics 
in each angular momentum sector /) depend on the choice of the original distribution of 
couplings A(Z) and of the total spin j. We remind that 



2j 



-n = WA(0) ■ A(0) = 53 |^-A(/) 



1=0 



2j + l 



A(l) = NA(l) 



(90) 
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and that Z(l) 

z(t) = M , A(i) = j2^' + iK-i) 2 ^' +l { 3 n J /')Hn (9i) 

A(OJ V=Q U J 1 ) 

In the trivial case where all the A's are equal we have 

A(Z) = A for all / A(0) = (2j + 1) N A and A(l) =0 if < I < 2j (92) 

Hence Z(l) = for all I > 0. We recover the GUE ensemble. 

We are interested in the situation where only a finite number Iq of interactions channels 
are activated, and where this number is much smaller than the total spin 

A(0 = if / > l , l < j (93) 

We first assume that all the A(Z) for I < Iq are of the same order A. Then A(0) is of order 
ANly/j. We thus rescale A — y A both with N (as before) and with the spin j, having in 
mind to have both a large environment (N —y oo) and a semiclassical spin (j —y oo). 

A(0 = (2j + l)^- 1 A(0 = (2j + 1)A(0 (94) 

and we are now interested in the limit of a large environment (N —y oo) and a semiclassical 
spin (j —y oo). 

— y oo , j < 1 large but finite or infinite , A(Z) of order 0(1) (95) 
We now have 

'° 

A(0) = t 2 = ^(21 + l)A(Z) independent of the spin j 

1=0 

In the figures presented in the Appendix [A] we plot for several choices of distributions 
of A's the resulting Z(l), and then discuss the results and prove some of the properties 
of the function Z(l). 



In fig. (16) we plot Z{1) as a function of I in the case where all the A(l) are equal 
for < I < 10, and zero otherwise. 



In fig. (17) we plot Z(l) as a function of I in the case where A(0) = 0, all the A(l) 
are equal for 1 < I < 10, and zero otherwise. 



In fig. (18) only a single and odd / contributes. 



In fig. (19) all odd I < Iq contribute. 



In fig. (20) all even I < Iq contribute. 



In fig. (21) we compare different cases, letting A(0) become large.. 
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• In fig. (22) we consider some random distributions of A(Z) for I < Zq. 

First we make the experimental observations: 

1. For < I < 2j, Z(l) is always in the interval ] — 1, 1[. 

2. For a fixed coupling distribution A, when the total spin j is large Z{V) takes a limit 
scaling form Y (which depends of course of the A's) 

Z(l)=Y(l/2j) when j ^ oo , l/j fixed (96) 

3. By normalisation Z(0) = 1. Otherwise Z(l) can be close to +1 (respectively to — 1) 
only in the limit j — > oo, when Z ~ 2j and when the A(/) are zero for all odd Z's 
(respectively all even Z's). 

3.5.2 Limit j — > oo, l/j fixed 

The existence of a limit distribution Y(x), x = l/2j when j — > oo is easily explained. We 



use Racah formula of the 6j-symbols to rewrite the formula 91 for A(l) as 



We use Stirling formula to take the limit 

j — > oo , x = — fixed (98) 
2] 

to obtain 

lo 

A{l)^J2 ( 2// + 1 ) F i' ( x ) (99) 

i'=0 



with Fi/(x) the polynomials 

(A;!) 2 (l'-k)\ 



F ^=i:7iSf^ x2k = ^(i+i',-!',!,^ (loo 



A:=0 

Hence the explicit polynomial form for the limit scaling function Y in eq. (96) 



£ A(0 (2Z' + l)F,,(a;) 

y(x) = ^ (ioi) 

E A(Z')(2Z' + 1) 
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3.5.3 Limit j — > oo, / <C j 

Apart from some very special cases, Z(l) is close to 1 only if Z is small. This case is needed 
for the study of decoherence. When I <C j we need only to keep the terms k = and 
k = 1 in the explicit form 97 for A(Z). We obtain 



A(Z) = A(0) - £ A(Z') (21' + 1) /'(/' + 1) + 0{{1/jY) (102) 

Hence 



1 



£ A(o (2z ; + 1) r(r + 1) 

Z(l) = l-l(l + l)--f- + ... , ^o = ^ (103) 

JU } £ A(Z')(2Z' + 1) 

This approximation is valid on the top of the curve F(x) near x = 0, i.e. provided that 

Z(Z + 1)< + (104) 

It is in particular valid when Z cx -y/J, which is the case to consider when studying the 
coherent states. 

Note that if all the A(Z)'s are of the same order A when Z < Iq and zero othervise, the 
coefficient D is of order 

Do~\l 2 (105) 

But note also that the numerator in D involves only the A(Z) for Z > 0. In particular, if 
the Z = coupling A(0) is much larger than the others A(Z) ~ A for < Z < Zo 

A(0)»Z 2 A =► ^o~^«Z 2 (106) 
3.5.4 The case of even Z's or odd Z's 



One sees from fig. (20) that when there are no couplings for odd Z's, i.e. when 

Z odd =^ A(Z) = (107) 

then for very large spin j — > oo the condition of slow decoherence Z(l) ~ 1 is satisfied for 
the Z's close to the maximal value Z = 2j. In fact in this case the scaling function Y(x) 
defined by eq. ( 96 ) is for x ~ 1 

Y(x) ~ 1 -cst.(l -x) (108) 

The reason why there is a very slow decoherence between opposite states such as \ j) and | — 
j) is of course that the coupling hamiltonian Hgs has an additionnal Z 2 parity symmetry, 
the inversion of spin, which protects the states with this symmetry from decoherence. 

On the contrary, one sees on fig. (18) and fig. (19) that when there are no couplings 
for even Z's, the condition of maximally fast decoherence Z{V) = — 1 is satisfied for the Z's 
close to the maximal value Z = 2j. 
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3.6 Decoherence parameters and norms of operators 
3.6.1 A(0) and the norm of H 

When all the A (7) are equal to the same A (GUE ensemble) all the matrix element of H 
are of the same order, vA, and the whole Hamlitonian if is a random M x M matrix in 
a GUE ensemble, with m = (2j + 1)N. The normalization is such that 



A(Z) = A tr(H 2 ) = (2j + l) 2 iV 2 A (109) 

In the general case where the A(Z) are different, we have 

tr(H 2 ) = (2j + l)iVA(O) with as before A (0) = — - J^(2/ + l)A(Z) (110) 

2j + 1 , 



so y A(0) is the typical size of an eigenvalue of H, that we call its norm \H\ 2 , the norm 
of an operator A being defined as r] 

II 112 tr ^^ V ) 

In our modified GUE ensembles, if we take for |m, a) a basis of the whole Hilbert 
space Hse — C^ +1 ' N , since 



(2j + 1)N 

m,n a,f. 



112) 



we see that our norm is 

\H\ 2 ~ x/(2j + l)iV x "average value" of \{m, a\H\n, (3)\ (113) 

3.6.2 The spectrum of H 



It is also easy to see from the solution 58 and 60 of the single resolvent H(x) that 



although the Hamiltonian H does not belong to a GUE ensemble but to the modified 
GU2,atE ensemble, in the large N limit its density of states (DOS) p(A) , i.e. the density 
distribution of the eigenvalues of the eigenvalues of H, is still given by a Wigner semicircle 
law, on the interval [-E, E] with 



p(X)ocyE 2 -X 2 with E = 2VA (114) 



2 This norm is the Hilbert-Schmidt norm divided by ^/dim(%), and is not the C* norm \\A\\ = 
su P|i/j) ^(v>|-0)^ which corresponds to the modulus of the largest eigenvalue of the operator. 
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3.6.3 Mean A av for the A(Z) 

It will be convenient to consider the mean value of the A(/) when averaged over all possible 
modes (Z,m)'s, A av . It is defined as 

^=(27TiF^ (2Z + 1) ^ (0 (U5) 

and in found to be nothing but 

Aav = 2^ A (°) = ( 116 ) 

The average value of Z(l) is 

(2i + l)^^— A (0) 

and in particular 

2j 



, , A(Z) /lif'LY 



H' = H- # (0) (118) 



if' being the purely 5 + £ interaction part of the Hamiltonian H, and being the 
purely external part £ of H. 



4 Decoherence and emergence of coherent states 

It is easy now to study the dynamics of decoherence and the emergence of the coher- 
ent states for spin as the semi-classical states robust against the interaction with the 
environment. 



4.1 Coherent states 
4.1.1 Pure coherent states: 

Coherent states are the pure states which minimize the uncertaincy relations for spin, i.e. 
the states with a maximally localized Wigner distribution. They read explicitely 

In) = V W 7 ^ - t cos{6/2) j + m sin(0/2y- m e- im< V) (H9) 

m t^- V (j+m)\(j~m)\ 

with (8, 0) the spherical coordinates of the unit vector n. Coherent states are formed by 
a coherent superposition of modes such that I ~ y/j <C j. Indeed for the single pure state 

K) = \j) (120) 
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the matrix density components are 

w ihm)_r w (i) w {i) __ I ((2j)!) 2 (21 + 1) 

w im - 5 m , w c _ s , , w c _ s - y (2j+z + 1)!(2j _ z)! ( 121 ) 

and for large j and small / Stirling formula gives 



The Wigner representation of the coherent state |n) is a Gaussian-like positive distribution 
with width 1/y/j centered at n on the unit sphere. 

4.1.2 Random pure states: 

At variance with coherent states, a random pure spin state \ip) is such that its density- 
matrix components are independent equally distributed random variables 

<& ~ 2JTT < 123 > 

and its Wigner representation is a random function on the sphere (analogous to a random 
polynomial with zeros obeying Wigner statistics with short distance cut-off l/y/J). 

4.1.3 Superpositions of coherent states: 

Quantum superpositions of coherent states correspond to more complicated functions. 
For instance it is well known that the "Schrodinger cat-like" state 

= j= (I + e z ) + | - e z » = { \ 3 ) + | _ j)) (124) 
corresponds to modes (/ even, m = 0) and (I = 2j,m = ±2j). 



ivW ! A h n n ((2j)!) 2 (2/ + 1) 1 /l + 2j 

^ CAT =2^,0 (i + (-i) )]] i2j + l + mj _ |+2*a- OW" 1 ) + y TT47 

(125) 

In general, the density matrix of a superposition of two coherent states will have large 
components describing the quantum correlations for 

W {lm) ^ for l _ j ( 126 ) 

In particular, if |^) is a superposition of two coherent states \ni) and |n 2 ) which are at 
distance 9 12 on the sphere S 2 , 

\ip) = Ci|ni) + c 2 |n 2 ) with n\ ■ Hi = cos6> 12 (127) 

the largest elements W^ lm ^ of the "off diagonal" part of the density matrix p fj = |n 2 )(ni| + 
\n\){n2\ occur for 

W^-0{1) for I ~ 2j sin(0i 2 /2) (128) 
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4.2 The time scales of the system 

From the above discussion, one sees that in our model the evolution of the spin is char- 
acterized by at least three time scales, which are such that 

r < n < r 2 (129) 



4.2.1 Dynamical timescale for the whole system tq = Td yn 

The first characteristic time is 

7-0= (A(0))~ 1/a = (|ff| 2 r 1 (130) 

It is the typical evolution time scale for a generic state of the whole system £ + S. Thus 
we denote it also Td yll . 



4.2.2 Decoherence timescale for spin states n = r^c 

The second characteristic time is 



n = r^V = vEh (131) 



\-Z m (\\W 



It is the time scale for the decay of the W^ m ^ coefficients of the density matrix for typical 
I ~ j. Hence it corresponds to the decoherence timescale for random pure states for 
the spin, with no particular spin properties. Hence we denote it also by t& cc . 

For t ^> ri, a pure random state \if>) has become a complete statistical mixture p oc 1. 
We have of course T\ > tq and in general T\ is of the same order than To. But we have 
seen that this decoherence time T\ can be much larger that r in the special case where 
A(0) ^> A(Z), I > 0, i.e; when the internal dynamics of £, given by = H^°\ is much 
faster than the dynamics of S induced by the coupling S + £ , and given by H$s = H' . 



4.2.3 Timescale for evolution of coherent states r 2 = r d 



iff 



The third characteristic time is the evolution time for a single coherent state. It is the 
decay time scale in sectors such that 



/~VJ«J (132) 
In this case the evolution of the density matrix is given by the regime Z(l) ~ 1 studied 



m 



3.5.3 Using eq. (103) this time r 2 is given by 

j 



T 2 = T 



Dr 



(133) 



with D given by eq. (103). Note that on general grounds 

— ^— ri < r 2 < jn where l = sup{/ : A(l) > 0} (134) 

As we shall see in sectionN |4.4 for t > r 2 the evolution of coherent states is decribed by 
a quantum diffusion process. Hence we denote it also by Tais- 
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4.2.4 The conditions for decoherence 



The dynamics of coherent states is much slower that the evolution of non-coherent states 



when the spin j is large and when l is of order 0(1), hence l < j. In fact using eq. (103) 



for D , and eq. (11) the ration of the two caracteristic times T\ and r 2 can easily be 



rewriiten in term of the ratio of operator norms 

\ 2 
111,) 

(135) 



r 2 < 'I 



I^la) (V' 



H' is the interaction Hamiltonian H' = ^ H^ l \ Of course the norm of the spin operator 

l>0 

S is simply \\S\\ 2 = + 1). The numerator is the squared norm of the commutator of 

the spin operators with H' 

(\\[s,H%) 2 = ^(ii[^,F']i 2 ) 2 = {2j l 1)N i2 tT [a s ^ H ']) 2 } ( 136 ) 

One sees explicitly that it is when the commutator of the spin with the interaction Hamil- 
tonian is small compared to the product of these two operators 

[S,H'] < SxH' (137) 

(here in the sense of the ||-|| 2 norm) that the single spin coherent states |ni) are much more 
robust than quantum superpositions of macroscopically distinct coherent states = 
Ci + c 2 |n 2 ) with # 12 yfj. It is in the regime 

n < t < r 2 (138) 

that coherent states behave as classical states (pointer states), and that one can use a 
semiclassical picture. 

In general, a pure state of the form = C\\ni) +c 2 |n 2 ) becomes a statistical mixture 
of \fii) and |n 2 ) after a time of order t\. But this decoherence timescale depends on the 
distance between the two states, i.e. of the angle 9 between the vectors n\ and f? 2 . In 
the semiclassical regime when the spin is large j ^> 1, and where this distance is small, 



but still large compared to the width of a coherent state, l/y/j <^ < 1, using eq. (128) 



eq. (103) and the scaling eq. (84) discussed in 3.4.4 , we have for typical decoherence 
timescale 

r de c.(0) ~ n/ sin(#/2) 2 such that n < r dcc .(0) < r 2 (139) 

Thus the closer two coherent states are in phase space, the longer it takes for decoherence 
to wash out quantum correlations between these two states. This is a well known ef- 



fect for harmonic oscillator or free particles discussed for instance in |Strunz et al., 2002 
Strunz and Haake, 2003 Zurek, 2003| . 



Beyond the regime discussed here, i.e. for t ^> r 2 , coherent states start to evolve, 
and in fact becomes statistical mixtures according to a quantum diffusion process that we 
discuss in section l4~4l 
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4.3 Illustrated examples of evolutions 
4.3.1 The cases considered 

We have a complete explicit solution of the evolution of the (/, m) components of the 
density matrix for the spin, starting from arbitrary initial conditions. The (I, m) com- 
ponents are nothing but the components of the decomposition of the Wigner transform 
of the density matrix in spherical harmonics. Therefore it is both tempting and easy to 
illustrate our results by explicitly plotting the time evolution of the Wigner transform for 
various initial states, and various choices of dynamics (given by the couplings A (7)), and 
various values of the spin j. 

In the following we take a large (but not tremendously large) value of the spin: j = 20. 
We chose as coupling distributions the simplest case 

A(0 = unless I = 1 (140) 

and express the time evolution in units such that tq = 1. This corresponds for the various 
time scales to 

7-0 = 71 = 1 , r 2 =j/2 = 10 (141) 

We represent the Wigner transform W(n) which is a real function over the unit sphere <S 2 
as a real function W(z) over the complex plane C using the stereographic projection 

n = (sin 9 cos (/>, sin 9 sin 0, cos 9) — > z = x + iy = re 1 ^ , r = 2 arctan(#/2) (142) 



4.3.2 Evolution of random states 

3 

We start with a random pure state oc ^2 a m \m), with the a m complex Gaussian 

m=-j 

independent random variables. In this case the initial W^ lrn ^s are (not independent) 
random variables of the same order (2j + 1) _1 , independently of I and m (except for 
tyr(o,o) _ (2j + l)" 1 / 2 ). For < t < n all the high angular frequency modes W^ lm) (t) with 
I ~ j undergo a fast decay, and for t > T\ only the low frequency modes with / ~ y/j are 
still there, which decay at a much slower pace. The evolution of the Wigner function is 



depicted in fig. (4.3.2) from t = to t = 8, that is approximately for < t < r 2 . One sees 
indeed the behaviour dicussed above. 

In other words, one starts from the pure quantum state \ip) which can be written as 
a quantum superposition of ~ 2j + 1 semiclassical coherent states. This state evolves 
and becomes for t > T\ a statistical mixture of these coherent states, with approximately 
equal probabilities of order 1/(2 j + 1) (with small fluctuations of relative order 1 / y/2j + 1 
which are responsible for the long-lived remanent wiggles. 



4.3.3 Evolution of a coherent state 



On fig. (10) we plot the evolution of a coherent state \n) 

\e z ) = \9 = 0,<P) = \j) . 



at the origin 



\9,4>), here the one centered 

(143) 
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Figure 9: Evolution of a random state for spin j = 20. The Wigner distribution W(n) is 
represented using stereographic projection of the complex plane (this is a .mp4 embedded 
file, use Adobe Reader® and click on the image to see the video). 



One sees immediately that its evolution in much slower and that it stays localized in 
phase space, as expected. We are working for finite j so that T2 is larger that r% but finite, 
which explain the fact that the Wigner transform widens with time. Its evolution will be 
discussed in more details in 14.41 
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Figure 10: Evolution of a coherent state for spin j = 20 (this is a .mp4 embedded file, 
use Adobe Reader® and click on the image to see the video). 



4.3.4 Evolution of superpositions of coherent states 

We now look at the evolution of quantum superpositions of coherent states. These kind 
of states are often dubbed "Schrodinger cat states" or "cat states" in the literature. 



2 state cats: On fig. (11) we start from a simple superposition of the two opposite 
(hence orthogonal) coherent states. 

l^> = ^(k/2,0) + |7r/2,7r>) (144) 

At time t = the Wigner transform shows the two peaks associated to the two coherent 
states and the strong interference fringes along the = ±7r/2 line (a great circle on the 
sphere), which contains the information about the superposition between the two coherent 
states. 

As time increases the two coherent peaks evolves slowly, while the interference fringes 
disappear very quickly with a characteristic time of order T\. We thus see explicitely the 
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Figure 11: Evolution of a "2-states cat" for j = 20 (this is a .mp4 embedded file, use 
Adobe Reader® and click on the image to see the video). 



decoherence induced by the couping of the spin with the environement, which transforms 
the quantum superposition of the two coherent states into a statistical mixture (here with 
equal probabilities) of these two states. 



3 state cats: On fig. ( 12 ) we start from a more complex superposition of three coherent 
states. 

|V>) = (|tt/2, 0) + |vr/2, 3vr/4) + [tt/2, 5vr/4» (145) 

At time t — the Wigner transform shows the tree peaks associated to the three coherent 
states and the interference fringes in between, which contains the information about the 
superposition between the three coherent states. The rightmost peak corresponds to the 
first coherent state \tt/2, 0), and the leftmost pair of peaks corresponds to the second and 
third coherent state |7r/2,37r/4) and |7r/2, 5tt/4). These two states are closer from each 
other than from the first one. 

As time increases the three coherent peaks evolves slowly, and the interference fringes 
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Figure 12: Evolution of a "3-states cat" for j = 20 (this is a .mp4 embedded file, use 
Adobe Reader® and click on the image to see the video). 



disappear quickly. One sees that the coherences between the close second and thirst states 
decay more slowly that the coherences between these states and the first one, as expected 
from eq. (139). 



4.4 Evolution of coherent states and quantum diffusion 



We now discuss in more detail the dynamics of a single coherent state at large times 
t > T\. We concentrate on the semiclassical limit where j ^> 1. Then the modes for a 
single coherent state are such that I ~ \fj <C j, the corresponding function Z(l) is always 
Z{1) ~ 1 and we may approximate the evolution function Ai^'\t) = M(t,Z(l)) by its 
universal form when 1 — Z(l) is small 



M {l \t) ~ W It 



4 D 1(1 + 1) 
3vr r j(j + 1) 



(146) 



with *f?(t) given by eq. (85) and 86 
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Using the asymptotics eq. (122) for the modes of a coherent state at large j, and 



starting at time t = from the coherent state \j) of eq. (143), the harmonics of the 
Wigner transform of the density matrix read explicitely 

W i. Ht)= 4±L exp (^)J t ±£i!^l±21) (U 7) 



There are two very different time regimes. 
4.4.1 Small time t < r 2 = r diff regime 

For the relevant modes I ~ y/j, the variable t' = t f°r ^ i s still very small, 



i! < 1, and in eq. g47|) ty(t') ~ 1. 

W®{t) = W«>(0) = exp (148) 

The coherent state does not evolve yet. The Wigner distribution in phase space W(n,t) 



given by 13 is a Gaussian distribution with width 



When t ~ r 2 the two terms W^(0) and ^f(t') are of the same order and the Wigner 
transform W(n,t) starts to widen in phase space. This is the regime we are seeing in 



fig. (10) since j is not very large. 



4.4.2 Large time t 3> r 2 = r^m regime 

When t ^> r 2 and I is small / ~ y/j the ^f(t') function dominates over the initial 
exp(— l 2 /2j) term. The coefficients of the Wigner distribution take the form 

V2JTT V 3tt r„ jU + l)J 
This gives a Wigner distribution W(n,t) in phase space (the unit sphere) with a width 



A e (t)~J— (150) 

The width of the distribution growths with the time t as y/i, thus the evolution of the spin 
in this semiclassical regime is a diffusion-like process. But as we shall see, the distribution 
in phase space is not exactly a Gaussian. 

4.4.3 Final time t 3> r cqu = jr 2 regime 

Finally, and as expected, there is a final equilibration time 

73 = JT 2 = r equ (151) 
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which corresponds to the equilibration time for the spin. When t 3> T3 = jr 2l only the 
I = mode survives, the quantum diffusion has completely averaged the spin over the 
phase space and we obtain a rotationally invariant complete statistical mixture of states 



W(n,t) = —==== (152) 
VM2j+l) 

It corresponds to the canonical distribution since the Hamiltonian for the spin is H$ = 0. 



We remind that from Sect. 2.2.2 in our normalizations for the Wigner distribution 



d 2 n W(n,t) = W-^— 
s 2 { J V 2j + 1 

4.5 Quantum diffusion and non-Markovianity 
4.5.1 The spin density profile 

We now look more closely at the diffusive regime r 2 <C t <C T3. Since the spin j is assumed 
to be large, a coherent state evolves into a mixed state whose Wigner distribution W(n, t) 
is still localized in phase space, with a width larger than the size of a quantum coherent 
state y/j, but much smaller than the size of phase space, of order unity. 
We perform the rescaling of the spin components 



n 



(u, Vl - v?) , u = y/2j + 1 z (153) 



so that the curvature of phase space becomes negligible, and phase space becomes the 
flat complex plane C with coordinate z = x + iy where z = (x,y). We perform a similar 
rescaling for the time 

t = t'r 2 = ^t' (154) 
and we normalize the Wigner distribution so that 

d 2 z W (z, t') = 1 (155) 



c 



In this large spin j and not too large time t regime, the phase space is similar to that 
of an harmonic oscillator and our results can be compared to those of decoherence mod- 
els for an harmonic oscillator, see in particular |Br aun et al., 2001 Strunz et al., 2002 



Strunz and Haake, 2003 



From 149 we find that in this limit the Wigner distribution is simply the Fourier 



transform of the universal decoherence function ^ 



W(z,t') = (2vr)- 2 J d 2 lm fe'l^T) e ife 

87rt' J_ n/2 V. t'cos(9)16J 
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Figure 13: The quantum phase space density profile W quan t- axn (z) (blue curve) versus the 
classical one Wa as ^ C ai(z) (red curve) in the large j limit. Normalization is such that the 
variance is the same(|2;| 2 ) = 1. There is a small but definite difference between the classical 
and the quantum cases. 



This result is valid of course as long as the width of the distribution is large but smaller 
that y/j (the radius of the rescaled sphere), i.e. for time scales such that 

72 < t < jt 2 = r 3 (157) 

In this quantum diffusive regime the Wigner distribution takes the universal scaling form 

W(z,t) = jW qiWDtum (r) , r = \z\/Vt (158) 



This universal function ^quantum i s plotted on fig. (13), together with the standart 
Gaussian which represents the density profile associated to a classical Markovian local 
diffusion process (Wiener process) in phase space. 

^1™! (z,t) = ^ eX P {-\zf / {it)) (159) 

4.5.2 Non-Markovianity 

It is interesting to notice that there is a small but definite difference between the two 
curves. Thus, once the decoherence caused by the coupling of the spin to the external 
system £ has absorbed the quantum correlations between the coherent states (and reduced 
in practice the quantum phase space to the semiclassical phase space, the sphere S2), the 
coupling between the spin S and £ still affects dynamics of the spin. This effective 
dynamics is a quantum diffusion process, not exactly a classical diffusion. Non trivial 
quantum correlations and quantum memory effects induced by the dynamics of £ and 
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the couplings are present at all time scales, resulting in the non trivial and non classical 
diffusion profile given by eq. ( 156[ ). 



Indeed, the self-similar profile W(z,t) given by eq. (156) can never be the solution of 



a homogeneous evolution equation local in time of the form 

— W(z,t) = M?W(z,t) , Ms differential operator in z (160) 
at 

(as would be the case if it could be written (at least in this limit t ^> r 2 ) as a quantum 
master equation of a Lindblad form). The evolution of the spin at times larger than 
the decoherence time T\ = t^ cc is described by a universal self similar but non-Markovian 



quantum diffusion process. As will be discussed more in section [6T2| such a non Markovian 
behavior is not so surprising. Non Markovian relaxation processes have been already 
observed and discussed in two level systems. 

This diffusion is universal in the following sense. It does not depend on the precise 
choice of the parameters A(Z) for the amplitude of the couplings in the different spin 
channels /, only from the fact that the couplings belong to a SU(2) xU(iV) GUE ensemble, 
and that they are important only for small /'s (in our case for a finite number of channels 
/ < /o). However the diffusion process cannot really depend on the internal dynamics of 
the environment £ (since this dynamics is given by the / = channel hence is GUE), nor 
of the initial state taken for S, since we take a random state iV'random) £ Hg. 



As we shall show in section |5.4[ a Markovian behavior emerges only in the limit where 
the dynamics for S becomes very fast, and when the initial state is chosen to be an energy 
eigenstate of for S. 

4.5.3 Spin relaxation 

Finally, as a simple application, let us compute the long time evolution of the spin com- 
ponent along one axis (here the axis z), assuming one starts from an initial spin state 
\ip s ) such that S z (0) = (ip s \S 3 \ip s ) ^ 0. From 

S z (t) = (S 3 )(t)=tr( Ps (t)S 3 ) (161) 

and averaging over the ifi n t ensemble we get simply the decoherence function for the I = 1 
mode 

S~ z (t) = M^(t) S z {0) (162) 

where (for large j) 

M^\t) t-^—^—] —-) (163) 



with ty(t) the function defined in section [3.4.4 So the magnetization relaxes to its equi 



librium value with a typical time scale r equ = r 3 , as expected. However, the relaxation is 
not exponential (as this would be the case in the quantum diffusion process were Marko- 
vian), but only algebraically as 

&(t) oc — , t oo (164) 

v 
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5 External dynamics for £ 



Up to now, and for simplicity, we have treated the internal dynamics of 8 exactly on the 
same footing than the coupling between S and 8. In particular the Hamiltonian 
(corresponding to the I = sector of H) was chosen to be a random Hamiltonian in a 
GUE ensemble, its distribution had the same U(N) invariance than for the / > sectors. 
No particular attention was given to the initial state IVw) f° r £, which was taken to be 
a random state. 



5.1 External Hamiltonian for 8 

Now the Hamiltonian Hs for the system 8, i.e. the I = component of the Hamil- 
tonian Hgs i s taken to be a fixed (not necessarily GUE) Hamiltonian with a regular 
normalized state density in the large N limit. We define the density of states p(E) and 
the normalized density of states v{E) by 

1 - 

u(E) = lim -p(E) , p(E) = E a ) (165) 

a=l 

and assume that v(E) is finite and regular. E a are the eigenvalues of the Hamiltonian 
H(°\ considered as a N x N matrix acting on Jig = C N . In the calculations we choose 
as a basis {\o), a = 1, • • • N} of T-Ls the eigenstates of H^°\ 

\a) = \E a ) , H<®\E a ) = E a \E a ) (166) 

For the S + 8 interaction (the I ^ sectors) we keep a random Hamiltonian H' in the 
SU(2)xU(iV) GUE ensemble, characterized by the variances 

A' = {A(Z), Z = l,---2j + l} (167) 



5.2 General form of the influence functional 

We can then repeat the calculation for the influence functional. We shall be interested in 
the evolution of the spin S, starting from a initial separable state of the form 

p(0) =p 5 (0)® \E a )(E a \ (168) 

i.e. an arbitrary spin state for S times a given pure energy eigenstate for 8. One can 
choose a more general fixed initial state \E) for 8, neither random nor a pure eigenstate 
of iff, but thanks to the Un invariance of the ensemble of the Hge it is in fact sufficient 
to consider the previously introduced case. 

The influence functional is now a tensor Ai (t, E) depending on t and E 

p S m(t) = tr £ (e-^(p*(0)®\a)(a\)e^) ru = M ru , st (t, E a )p s st (0) (169) 
3 Strictly speaking =1 S ®H £ 
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with H = + H' . We shall compute this functional in the large N planar limit. 
We better consider the double resolvent 



G?Ct(%,y) = (H 

and the single resolvent 



x-H 



y-H 



\ua) 



1 



\sP) 



[ x-H< 

Obviously, when there is no coupling between S and £, this resolvent is simply 



ra 



x 



1 



x — E n 



(170) 



(171) 



(172) 



When there is a coupling, i.e. when the A(Z)'s are non zero, thanks to the SU(2) 
invariance of the distribution ensemble for the H' and the fact that the \a) are eigenstates 
of H(°) we can show that 7i"f(x) is still of this form 

H?f(x)=5 r>s 8 a ,pC a (x) (173) 

In the large N limit, C a (x) satisfies the recurrence equation (which generalizes the recur- 
rence equation 59 for T-L(x) ) 

C a (x) 

where 



D 



N 



C a (x) \ 



X 



(174) 



2.1 



5 = £^A( 



(175) 



^ 2 j + 1 

i=i J 

Thanks to the SU(2) invariance, we can easily show that the double resolvent Q is 
diagonal. Its double Wigner transform decomposition takes the form 



?<x,P 



;i76) 



Q rn (x,y) is given by a planar recursion equation similar to eq. (69) for Q ^ (x,y). Its 
solution is 



( x ,y) = 5 a ,pC Q (x)C a (y) + 



D{1) C a {x)C a {y)Cp{x)Cp{y) 



N 



1 - D(l) [ J2 C 7 (x)C 7 (y) 

V7=l 



(177) 



where the D(l) are defined, as in eq. (66) (same SU(2) structure) by 



2j 



/3(/ 1 ) = ^A(/ / )(2/ / + l)(-l) 2j+r 



i'=i 



j 3 h 



(178) 
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Note however that in eq. (178) the sum over I excludes the I = case, contrarily to the 
sum in eq. (66) which defines the A(/)'s. 

Since we are only interested in taking the trace over the final states \/3) we simply 
have to consider 



N 



C a {x)C a (y) 



7=1 



(179) 



It is of course natural in the large N limit to rexpress the sum over states of £ as a 
continuum integral over the spectrum of H^°\ 

N r ~ ~ 

-> N / dEv{E) , C a {x) -> C{x,E a ) , Fgfoj/) F (l) (x,y,E a ) (180) 



eq. (174) becomes 

C(x,E) 



with 



x - E - A'C(x) 



2j 



A' = iV/3 = ^^±iA(/) 



^ 2 j + 1 



dE v(E)C{x,E) 



A(0 = iVA(Z) 



1811 



182) 



while eq. (179) becomes 

F(l)(x,y,E) 

where 



C(x,E)C(y,E) 



1 - A'(Z) / dE' z/(£') C(x, E')C(j/, £') 



2j 



A'(/) = ND(l) = ^A(r)(2/' + l)(-l) 2j+ '' 



i'=i 



183) 



184) 



The notation A' and A' (I) (with a tilde) in the definitions of eq. (182) and eq. (184j) is 
here to recall that there is no I = contribution in the sum over the Vs, contrary to the 
definition for A(Z) and A(0) given by eq. (661) and eq. (67). We have obviously 



A'(0) = A' 



(185) 



As in the previous situation, the evolution of the reduced density matrix factorizes 
into each (l,m) sector, and the influence functional becomes a single function of the time 
t, of the angular momentum /, and now of the initial state energy E. Thanks to the SU(2) 
invariance it is still independent of m. It is given by the integral 



M (l \t,E) 



dx f dy _ it 



2m T 2ivr 



Fm(x,y,E) 



(186) 
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This can be rewritten in a simpler form by the change of variables (and its inverse) 

w = W(x) = x- A'C(x) , x = X{w) (187) 



Indeed eq. (181) becomes 
C{x) 



dE 



u(E) 



C(x,E) 



(188) 



w — E w — E 

so that C(X(w)) is the Hilbert transform of v(N), i.e. the resolvent of H^. Then in 
eq. (183) we rewrite the integral 

dE' u(E')C(x 1 ,E')C(x 2 ,E') = J dE'u(E') 
C(W(x l ))-C(W(x 2 )) 



1 



W(xt) - W(x 2 ) 



1 
A' 



W(x 1 )-E'W(x 2 )-E' 

X\ — X2 



W(xx) - W(x 2 ) 



+ 1 



(189) 



and after some algebra the evolution kernel Ai^(t,E) is written in a form similar to 
eq. (77) as an integral representation involving the variances A(Z)'s for the interaction 
Hamiltonians H$e — s ^2 i H' yl,m ' > and the Hamiltonian for E, He = through the 

function W(x) (related to the inverse of the resolvent for Hs)-, an d a parameter Z'(l) 



M {l) {t,E) 



dx\ f dx 



-it(xi—X2) 



where 



2m J 2m {W{x x ) - E)(W(x 2 ) - E) (1 _ Z , {1)) + Z , {1) ^ 

A'(i) 



-X2 



Z'{1)) 



A'(0) 



)-W(x 2 ) 

(190) 
(191) 



with A' (I) defined through eq. (178) and eq. (183). 



5.3 Application to the Wigner ensemble 
5.3.1 General form of the solution 

As a simple and illustrative ensemble, let us treat the case where the density spectrum of 
Hp is the Wigner semi-circle distribution with width 2E . 

v{E) = ^rjEl-E* (192) 

This corresponds to take for H^ a sample in a GUE distribution, as in the previous 
section, with the identification 

r — N 

E = 2^A(0) where A(0) = ^yA(O) = A av (193) 
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But now we are able to study the dependence of the final state on the choice of initial 
state \E). 

To compare our calculations with the results of the previous section, it is convenient 
to perform the same rescalings. The first time scale r is now 



To 



A(0) 



A(0) = A av + A'(0) = ^ + f:^A(/) 



^ 2 j + 

1=1 J 



(194) 



The Z(V) parameters of eq. (76) (not to be confused with the Z'(l) of eq. ( 191[ )) and the 
average Z av of eq. (117) are 



Z(l) 



A av + A'(Z) 



A s 



As, 



A av + A'(0) 
After performing the rescaling 

t = t/r , 

the evolution kernel MS l,m \t, E) takes the form 

M {l ' m \t,E) = M(i,E,Z(l),Z, 



A(0) A av + A'(0) 
E = Et 



(195) 



(196) 



(197) 



The function M(t, E, Z, Z av ) is given by the integral representation, which generalizes 

e -it(X 1 -X 2 ) 



eq. m 



M(t,E,Z(l),Z a 



dX 1 f dX : 



Z eM H x H 2 



2ivr J 2ivr (W 1 - E)(W 2 - E) 1 - Z{$)R X H 2 



(198) 



where as above H\ i are given by the function H(X) (the resolvent for the Wigner ensem- 
ble) 



1 



1 



H h2 = H(X h2 ) , H(X) = -(X- VX 2 -4) i.e. X = H + — (199) 

2 n 

while Wip are given by the function W(X) defined from H(X) and Z av through 

fW(X)\ 



W li2 = W{X lfi 



'Z^H(X) = H 



i.e. W{X) = Z av H(X) + 



H(X) 
(200) 

As in the previous section, H(X) is analytic in C\[— 2,2], with the cut along the 
interval X G l\ = [—2, 2] corresponding to the spectrum of the total Hamiltonian H = 
Hg + Hse- Away from this cut, one has X G C\Ji \H(X)\ > 1. In addition 

H(W) is analytic in C\[— 2\/Z av , 2\/ Z av ], with the cut along the interval W G I 2 = 
[— 2y/ Z av , 2sJ Z av ] , included in the previous cut since Z av < 1. Away from this cut, one 
has W G C\/2 =^ \H\ > ly/Z av . Remember also that E is the (rescaled) energy of the 
initial state \E) G Tig, chosen to be an eigenstate of H?, hence \E\ < 2^Z av . Therefore 
the integral form for M(i, E, Z(l), Z av ) given in eq. (198) is defined for a contour in X\ 
and X 2 encircling the cut I\ = [—2, 2] and the poles W (X) = E, or equivalently for a 
contour in H\ and H 2 such that \H\ > l/y/Z av . 
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5.3.2 Discussion of the solution 



Altough it depends on more parameters, the function M(i, E, Z, Z av ) has a large t behavior 
similar to the behavior of M(t, Z) depicted in fig. ([6]). As long as Z < 1 it decays when 
t — > oo (fast decoherence) . 

In particular, if we choose Z av = and E = we recover the function M(t, z) given 
by eq. (77), depicted on fig. ^ and studied above. 

M(t, E = 0, Z(l), Z av = 0) = M(t, Z(l)) (201) 




Figure 14: The decoherence function M(t, E, Z, Z av ) plotted as a function of t and of Z 
for E = and Z = Z av . Its behavior as a function of time t and of Z av is similar to that 
of the function M(t, Z). 



For typical Vs (i.e. I ~ j) Z(l) ~ Z av . If we choose 



z(0 



Z av and £/ = 



we get 



M(t,0,Z., 



iV(t,Z a 



ciX e 



-ax 



(202) 



(203) 



2ivr W(X) 

The function iV(t, Z) is easily calculated, for instance from its small t and small Z series 



expansion. The function M(t, 0, Z, Z) is depicted in fig. (14). Although not exactly the 



same as the decoherence function M(t, Z) obtained in the previous section (see eq. (74)) 
and depicted on fig. ([6]), one sees that it is very similar. In particular one sees that when 
Z av is close to 1, the decoherence becomes very slow. 
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5.3.3 Large t asymtotics 

This is a general feature: for fixed Z av < 1 and fixed E, when I is small (i.e. I -C j), Z(l) ~ 
1 and the decoherence becomes always very slow. This can be shown by reexpressing the 



integral representation 198 for the function M(i, E, Z(l), Z av ) by a double contour integral 
over Hi and Hi- The contours can be choosen to be the unit circle = e 1,2 , and at 
large time t the integral can be evaluated by the saddle point method. The relevant 
saddle points are the extrema of the functions X(H) = H + 1/H, situated at Hi and 
H 2 = ±1. We shall not detail the full calculation here, but one can show that the final 
large i behavior of M(i, E, Z(l), Z av ) is of a similar form than for the function M(t, Z(l) 



given by eq. ( 82 ) . 



M{t, E, Z{1), Z w ) = t~ 3 (A — B sin(4t)) (l + 0(t _1 )) (204) 

The coefficients A = A(E, Z(l), Z av ) and B = B(E, Z(l), Z av ) depend in a complicated 
way of the spin / through Z(l), but now also of the energy E and of Z av . 

Let us just give the form of these coefficients in the limit of small angular momentum I. 
This is the limit where the decoherence becomes slow and where coherent states emerges 
as semiclassical states. In this case we have seen that Z(l) is very close to 1, and we find 
that in this limit 

~ 7T (1 + Z av )2 _ £2 (1 _ Z{ l )f \ ) 

B ~ 0(1) (206) 

One sees in particular that in this regime the function is a simple scaling function of the 
variable t(l — Z{1)). 

5.3.4 Small Z(l) scaling 



One can study the regime (similar to the regime of 3.4.4) where 



Z(l) -)• 1 , t(l - Z(l)) = t' = 0{l) (207) 



In this limit the dominant contribution in the integral 198 is given by H 2 — \/{Z(l)Hi). 
Noting Hi = e~ w and integrating over H 2 we find the scaling function 

M 8ca n n g(t, E, Z(l), Z ay ) — j —— — Z "Z! ,2, — ~ e_z ' " ' 



'o t + z m ) cose- e) 2 + (l- z^y (sine) 2 

(208) 

For large t' this integral is dominated by the contribution of the endpoints 9 = 0, tt, and 
this function scales indeed as 



M r ~ 1(1+^+^(1^ (2Q9) 

M scalmg _ ^ (1 + Zav)2 _ j£;2 f ,3 



However, when 1 — Z av is small there is an interesting regime where this integral is 
dominated by the point 8 C where the denominator is small, namely 

(l + Z av )cos# c -E = (210) 
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The integral can be approximated by a Lorentzian integral and one finds that the function 
■A^scaUng scales as 

(211) 



M 



scaling 



exp -t' \JA-E 



The crossover between the exponential decay 211 (valid for t' < t' CTOSS ) and the algebraic 



decay 209 (valid for t' > t' CIOSS ) occurs for 

4 oss ~ io g (i/(i - z„)) = login/ro) 



(212) 



5.4 The limit of fast £ dynamics (n ^> tq) 
5.4.1 Fast dynamics 

These results can be used to study the regime where the dynamics of £ is much faster 
than the dynamics on S induced by the coupling S + £ . In the situation considered here 
(Wigner ensemble for H?), this occurs if the width of the spectrum of H?, namely 2En , 



is much larger than the width of the spectrum of H$£, which is of the order of \ A'(0) 



En 



2\/A av > VA'(O) 



(213) 



This situation is depicted on fig. (15). 




Figure 15: The case of fast £ dynamics. The energy spectrum density u(E) of is 
depicted in blue. E is the energy of the initial energy eigenstate \E) of £ that we couple 
with the spin S at time t = 0. The energy perturbation of this state and its mixing with 
other energy eigenstates \E') of £ induced by the coupling Hamiltonian H$£ is represented 
in red. 



In this case we use the results of Section 4.2). The time scale tq for the internal 



dynamics of the whole system is from eq. (194) 



TO 



A(0) 



A a 



2 



(214) 
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The decoherence time scale t\ for a typical spin state is from eq. (131) 



r 



1 



r 



A'(0) 



hence it is much larger that tq. 



We remind that from eq. (194) A'(0) is defined by 

A'(0) = X>1 + 1) A(Q 
z=i 



(215) 
(216) 

(217) 



where the A(l) are of order 0(1). The factor Z(l) behaves for small I according to eq. ( 103 ) 

Z(l) 

with now 



1 1(1 + I) Do 
J(J + 1) 4 + 



(218) 



-■ '0 

A) - — ^A(Z)Z(Z + 1)(2/ + 1) < 1 



(219) 



av j/ =1 



We have seen that the decoherence time scale for a coherent state is from eq. (133) 

3 



T~2 = T 



Dr 



(220) 



Therefore when the spin j is large it is always much larger than T\ since 

T2__ jEfii(2/ + l)A(/) j 



~ — ^> 1 

• n 



^ E^i'(i + l)(2i + l)A(0 
5.4.2 Evolution of coherent states 

In the limit of large Eq, from eq. (213) and eq. ( |195 ) we have in any case 

Z w ~ 1 and Z(Z) ~ 1 



(221) 



(222) 



and we can apply the results of the section 5.3.4. Let us consider as an initial state |^) 
of the system S + £ the single spin coherent state \e z ) times the Hs eigenstate \E). 

\9(t = 0)) = |4) ® |^) (223) 

We study the evolution of the reduced spin matrix density for times t ^> r , as was done 
in Section 4.4 for the case of a spin coherent state \e z ) times a random £ state | 



. The 

harmonics of the Wigner transform of the density matrix are now, using eq. (211 ) for the 
scaling of the evolution functional M(t, E, Z(l), Z av ) 

21 + 1 f l 2 \ ( . 1(1 + 1) D 



W®(t) 



V/2JTT 



cxp 



2j 



cxp 



-t 



3(3 + 1) 4 



El 



E 2 



(224) 



to be compared with eq. (147). This gives a Wigner distribution W(n,t) on the unit 
sphere with width proportional to \/i when t T2, hence the evolution of the spin is 
again a quantum diffusion process. 
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5.4.3 Markovian limit for quantum diffusion in phase space 



For fast £ dynamics, initial state \E) and large spin j > 1, the spin diffusion at t 3> r 2 
becomes Markovian, with no memory effects. Indeed, at not too large times, whan we can 
approximate the phase space 5 2 by the tangent complex plane C, the Wigner distribution 
W(u,t), which is the Fourier transform of the decoherence function 



exp 



E 2 



is exactly a Gaussian distribution of the form 159, given by eq. (156) 

1 



W{u,t;D) 



AnDt 



exp(-|u|7(4Dt)) , 



(225) 



instead of the non-classical distribution function Wquantum of eq. (156). This corresponds 
to a classical diffusion process (Wiener process) on the unit sphere^] S2 with a diffusion 
coefficient D which depends on the energy E of the initial state through 



D 



El 



E 2 



(226) 



This Markovian regime where quantum diffusion is classical is expected to hold only when 
the exponential scaling 211 of M sca ii ng is valid. Since it holds only for t(l — Z(l)) < t' cr:oss 



given by 212 , we can show that non Markovian deviations to the distribution function 
W(Z,t) given by 



225 



occurs when (this is a very rough estimate) 



\u 



> AD log(n/r ) 



(227) 



This means that departure to Markovian behavior for finite T\ j tq are more easily observ- 
able in the large distance/small time diffusive regime. 

5.4.4 A Golden Rule formula for the diffusion coefficient 

It is interesting to note that the diffusion coefficient D can be rewritten in term of the 



normalized density of states v(E) for the system 8, given by eq. (192) 

u(E) 2 
and of the norm of the commutator 



nE% v ~~ ^ 



C = i[S,H se ] 



(228) 



(229) 



between the spin operator S and the coupling Hamiltonian Hse- Indeed, using eq. (219) 

c 

1 



for Do, eq. (193) for Eq, and the eqs. (133 137) of sect. 4.2.3 we have 

1 



D = 2nu(E) 



N(2j + i; 



rtr [C 



2ttis(E) 



2j + 1 



CI 



(230) 



4 Important: we do not perform the rescaling 153 by a factor y/2j + 1 on the spin here. 
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Using eq. ( |165 ) we can rewrite the diffusion coefficient D in term of the full density of 
states p(E) = Nu(E) for the Hamiltonian Hg, acting on the full system, and of the 
average squared norm of a matrix element of the commutator C 



E 



;#|C7|*') 



N(2j + 1) 



CI 



(N(2j + 1)Y ^ 



Y,\mcm (23i) 



where {|^)} is a basis of H$g = C N ^ +1 \ Indeed we rewrite eq. (230) as the general form 
for the spin diffusion coefficient D 



D = 2np(E) E \(^\C\^' 



(232) 



This form for the diffusion coefficient for the dynamics of coherent states, as a function 
of the energy E of the initial state for the environment, is similar to the Fermi Golden 
Rule for the transition rate of a state if) of a quantum system into a continuum of states 
\&' (with density of final states pfi na i), induced by a perturbation Hamiltonian H pert . 



2np & 



nal 



ffpertl* 7 )! 5 



(233) 



This is not so surprising since we are in a regime where effective dynamics is Markovian, 
and it is known that there is a relation between the Markov approximation and the 
Fermi golden rule for simple systems coupled to a bath (see |Alicki, 1977| and references 
therein). The form 232 for the diffusion coefficient can therefore be viewed as a fluctuation 
dissipation relation for the quantum dynamics of the spin in the regime t ^> ti. What is 
interesting is that it involves the matrix elements of the commutator C between the spin 
operator S and the interaction Hamiltonian H$g- 



We therefore expect that the relation 232 for D is general and independent of the the 



explicit form of the density of states p(E) for the external Hamiltonian Hg. It should be 



easy to study the limit t\ 3> T2 using the general form 190 for the evolution functional 
M(t, E, Z(l), Z av ) instead of the explicit form 198 for a semicircle state density law. In 



particular, we remark that the derivation of the exponential scaling 211 starting from the 



general form 208 of the scaling function M sca i ing at small I uses a Lorentzian approximation 



very similar to the one used in the textbooks derivation of the Fermi Golden Rule. 



5.4.5 Initial state \if>s) dependence and randomized Markovian processes 

It remains to understand why, even in the limit of fast dynamics T\ 3> To, if for £ as 
initial state l^w) = 1^(0)) m Hg one starts from a random state |V> ra ndom), i n the quan- 
tum diffusion regime the distribution profile in phase space W(u,t) for the spin is not a 
Gaussian, but is given by the non-Gaussian quantum distribution Wquantum of eq. (156) 



given by the Fourier transform of the scaling function \I> obtained through eq. ( 85 ) 



This can be understood easily. First let us come back to the case where Hg is a GUE 
random hamiltonian, so that its normalized density of states v{E) is given by the Wigner 
semi-circle law 12281 

HE) 



kE% 



El 



E 2 



(234) 



50 



In the fast £ dynamics limit T\/tq — > oo and large spin limit j — > oo, if one starts form 
a separable initial state product of a coherent spin state \n) and of a random £ state 
| ^random) , in the quantum diffusive regime r 2 <C t <C jr 2 the probability distribution for 



the spin 156 can be written as 



W(u,t)= dEu{E)W{u,t;D(E)) 



(235) 



with v(E) given by 234, W(u,t; D) the Gaussian density profile 225 and D(E) the spin 
diffusion coefficient 226 when starting from an initial energy eigenstate \E). 



W{u,t;D(E) 



AnD(E)t 



exp 



AD(E) t 



D(E) 



£ 2 



4j(j + 1) 



(236) 



Eq. 235 has a simple probabilistic interpretation. The spin diffusion process is simply 
a statistical superposition, with probability distribution u(E), of the Wiener processes 
with the different diffusion constants D(E). In other word, starting from a random state 
for £, to estimate the probability (density) to observe the spin at position u in phase 
space at time t, for all practical purpose we may assume that: 

1 - at time we choose at random an energy eigenstate \E) for £ (therefore with proba- 

bility law given by the density of states u(E)), 

2 - then we let the spin diffuse according to a Wiener process with diffusion constant 

D(E). 

This kind of non-Markovian process is build by the operation of randomization of a fam- 
ily of Markovian processes (see the mathematical litterature, for instance |Feller, 1971| ). 
Here the family of Wiener processes depending on the parameter E through the diffusion 
constant D(E) is randomized with the probability distribution v(E). 

To understand the physical origin of this ramdomization, we remark that a random 
state | ^random) of £ is nothing but a quantum superposition of energy eigenstates of Hg 
with random complex coefficients ce, choosen independently according to the Gaussian 
normal distribution. 

^random) = Y, CE \ E ^ ( 23? ) 
E 

Hence the initial step 1 amounts to project the initial state ^random) onto a random energy 
eigenstate \E), i.e. to perform an ideal (von Neumann) measurement of the energy of the 
external system £. 

This is to be expected. Indeed, we are in the regime where the interaction Hamiltonian 
Hss is a small perturbation compared to the internal Hg Hamiltonian for £. But we are 
also in the large spin regime where j ^> 1, so the spin S can be considered as a large 
quantum system, with many states, weakly coupled to £. In the same way than spin 



decoherence occurs by the time scales T\ with (from 130 and 131 ) 



TO 



\H. 



sej 



\H. 



se 



\H, 



< 1 



(238) 
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with H = Hse + Hg, we expect that the decoherence between (not too close on the 
spectrum E line) energy eigenstates \E) of Hg occurs at a time scale r e given by 



To 



\H\\ ■ \\Hg 



\[Hs£,H £ ] 



< 1 



(239) 



In other words, the interaction between S and £ induces both decoherence for the spin 
and decoherence in energy for the external system. As a consequence, at the time scale T2, 
where quantum diffusion starts, the initial coherent state \n) for the spin can be considered 
as a classical spin state n, while the initial random state | ^random) for £ can be considered 
as a classical equidistributed statistical mixture of energy states E. 

This suggests the general form of the distribution profile W(n, t) as a function of 
the £ initial state |VWt)- We consider a general Hg with a continuous density of states 
p(E) = Ni>(E), and we normalise the energy eig enstates |E) = (l/y/p(E))\E) so that 



(E|E') = 5(E - E') 



(240) 



If we start from a general initial state \ip) for £ which has a smooth decomposition in 
energy eigenstates, that is to say 



(E|-?/>) = ip{E) is a regular function of E 



(241) 



at least when averaged over small energy intervals AE ~ 1/ T 2> then the probability 
distribution for the spin at large time t 3> r^g should be given by 



W(u,t;4>) 



dE\(E 



AttD(E) t 



exp 



u 



AD{E)t 



with the effective diffusion constant D(E) given by eq. (232) 



D(E) = 2vrp(£)E 



E = E 



tr [(i[H S£ ,S}y 
(N(2j + l)) 2 



(242) 



(243) 



6 Discussion and conclusions 



6.1 Summary 

We have proposed and solved a model of a general quantum spin j interacting with a 
large external system (with iV states). The interaction Hamiltonian H mt is described by a 
new random matrix model, relying on a GU2x7vE ensemble, which takes into account the 
interactions in all possible angular momentum channels. The only constraint is that the 
probability distribution for H is Gaussian and invariant under the spin symmetry SU(2) 
and the external U(iV) group. The dynamics if the model is solved exactly in the large 
N limit, in the case where the internal dynamics of the spin is trivial. We can thus study 
in full generality the decoherence of spin states in various dynamical regimes, and as a 
function of the initial condition and of the dynamics of the external system. 
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In the semiclassical j — » oo limit, we are able to characterize quantitatively the interac- 
tion Hamiltonians which are such that coherent spin states have a much longer coherence 
time than any other states, and emerge effectively as the semiclassical "pointer states" of a 
classical spin. The large time dynamics of these coherent states is found to be a quantum 
diffusion process, in general non-Markovian. 

When the dynamics of the external system is much faster than the dynamics induced 
by its coupling with the spin, and when the initial state of the external system is an energy 
eigenstate, the quantum spin diffusion is found to be Markovian. The diffusion coefficient 
takes a golden rule form, but involving matrix elements of the commutator between the 
spin and the interaction hamiltonian. When the initial external state is a general state, 
the quantum spin diffusion is described by a randomization of Markov processes. We 
argue that this is explained by the decoherence of energy states (for the external system) 
induced by the coupling with the spin. 



6.2 Comparison with previous studies of decoherence in spin sys- 
tems 

We finally discuss the relation of the present work with some previous studies of open 
quantum spin systems. 

In the pionnering work |Takahashi and Shibata, 1975 , the relaxation of a large spin j 



in an external magnetic field and coupled to a heat bath is already studied. A coherent 
state representation for the spin states different from ours is used. The Born approxi- 
mation and some methods of |Agarwal, 1971| are used. Therefore their results are only 
valid for large time. A relaxation towards a classical distribution is obtained and the 
corresponding diffusion coefficients are calculated. In |Shibata and Saito, 1975 the spin 



relaxation in the specific case of the j = 1/2 spin is studied within the Markovian ap- 
proximation. In both papers the spin has a non trivial dynamics since it is coupled to an 
external field, this in not taken into account in our study. 

In |Mello et al., 1988| the relaxation of a j = 1/2 spin coupled to an external large 
system (bath) is studied. The spin has no internal dynamics (H sp [ n = 0), the bath has 
its own dynamics characterized by some Hb with its density of states. The interaction 
is of the simple form H int = a x ® V ext , where V ext is a random operator for the external 
system, chosen in a GOE ensemble local in energy E (for instance the matrix elements 
of V e ^t are non zero only between energy eigenstates such that the energy difference is 
finite \E\ — E%\ < AE, thus it is a finite width band GOE). They use the RMT methods 
which shall be used in most subsequents studies based on RMT, including the present 
one. With the initial condition considered they are able to study the relaxation of the 
spin at all time, and find an exponential relaxation (with already a Golden rule form for 
the relaxation time). 

In |Lebowitz and Pastur, 2004] and Lebowitz et al., 2 007 the study of relaxation for 



this spin j = 1/2 model with a band GOE H[ nt has been extended to the non trival 
case where the spin has its own dynamics i^ spin oc o z . Special attention is paid to the 
case where the external bath is formed of large number of independent large subsystems, 
in some particular limits where an effective temperature can be associated to the large 
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bath. Explicit expressions are obtained for the evolution of the spin 1/2 density matrix, 
but in practice the relaxation regime towards an equilibrium distribution is studied, in 
particular in the so called van Hove limit where the coupling spin-bath is small, which 
corresponds to our limit r dyn <C r dec . Thus the study of [j Lebowitz and Pastur, 2004 



Lebowitz et al., 2007| is more limited than our work, since only the spin j = 1/2 case is 



considered and the GOE ensemble for H\ nt is not SU(2) invariant, but also more general, 
since H mt is a band matrix in the energy states for the bath, and, most importantly, the 
spin is coupled to an external magnetic field, so that relaxation and thermalization effects 
can be studied. 

Another study involving RMT techniques is |Lutz and Weidenmiiller, 1999 . General 



system+bath couplings with band random matrix Hamiltonians are studied, both for a 
quantum harmonic oscillator and for a two level system. In the case studied, Markovian 
equations can be derived. 

In |Esposito and Gaspard, 2 003 an extensive study of the relaxation of a two level 
system (spin 1/2) coupled to a bath is performed. The coupling is given by a H mt = 
o'x <8> H>ath Hamiltonian with Vbath a GUE random matrix, ifbath also given by a random 
matrix, and H spm oc a z . Most of the study is by numerical simulations. The finite iV 
effects, where N is the size (number of states) of the bath, the statistics of the eigenvalues, 
the crossover between the Poisson and Wigner behavior for the spacing of the eigenvalues 
are very thoroughly studied. Explicit checks of the self averaging property for different 
realization of the random Hamiltonian are done. The relaxation and equilibrium effects 
are considered. Interesting analytical results are obtained in the regime of weak coupling 
(which corresponds again to the regime Td cc 3> 7d yn in the present paper). In particular 
the fact that the relaxation can be non-Markovian is observed and discussed. The fact 
that in the regime of strong coupling the decay of the magnetization (evolution of (cr z )) 
with time can be non-exponential but algebraic with some oscillations is also observed. 

The relaxation of a two level system (and systems with a small but > 2 number of 
levels) coupled by RMT Hamiltonians to a bath have been also studied numerically and by 
the so called HAM method in |P. Borowski et al., 2003 Gemmer and Michel, 2006] . The 



two level system has been studied more extensively and precisely by the TCL method in 
|Breuer et al., 20 06 j. Here also, it is found that depending on the initial states for the 
system and the bath, the relaxation dynamics can be non-Markovian. 

Finally, let us quote some papers which discuss more specifically the problem of de- 
coherence in spin systems. There are of course many studies of the dynamics of several 
coupled spins 1/2, and of their equilibration dynamics, see for instance |Saito et al., 1996 



In a recent numerical study |Yuan et al., 2009] it is confirmed for instance that decoher- 
ence is a much faster process than the thermal relaxation. 

In |Braun et al., 200 1 |Str unz et al., 200"2] |Strunz and Haake, 2 003 decoherence is 
studied for various systems (mostly a particle or an oscillator) coupled to a bath through 
a simple H int = Q sys t <8> -Bbath Hamiltonian (sometimes two), in the so-called interaction 
dominated regime, which should correspond in the present paper to the regime r syst <C r<i ec . 



Of particular interest for the present work is the section VII of |St runz et al., 2002 , where 



the system considered is a large spin j, studied by coherent states techniques. However 
the difference is that the coupling agent is just Q sys t = cr x , while the bath Hamiltonian 
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Sbath is a sum of in dependant Hj's for small independant sub-baths (so that random 
limit theorems may be applied). Thus it is not clear how to compare this model to our 
model. Some small time expansion are used and with these approximation the decay of 
the off-diagonal elements of the density matrix is obtained (decoherence) but with a time 
dependence which is very different from the exact results that we have obtained. 



6.3 Generalizations and open questions 

Many interesting questions remain to be addressed for this kind of models of a spin coupled 
to a large external system. Firstly, we have seen that the non-Markovian dynamics that 
we obtain for large spin j has been already observed in some two level systems. In order to 
better understand this effect, and how the non-Mar kovianity is related to the dynamics 
of the coupling and to the initial states, one must generalize the random interaction 
Hamiltonians that we have considered to the more general case of an gaussian ensemble 
with a law depending on the initial and final energy states for the external system, as in 
|Mello et al., 1988 Lebowitz and Pastur, 2004] . It should not be difficult to extend the 



study of section [5] to this more general case, but this will be discussed in a further work. 

In order to compare more precisely our results with those of previous studies, one 
would like to study simpler coupling Hamiltonians of the form H int = Q sys t <8> -Bbath with 
for instance Q sys t = S x . However, our model is solvable and its solution takes a simple 
form precisely because the random interaction Hamiltonian belongs to a ensemble H iat 
which is SU(2) invariant. This is clearly not the case for the simple .Hint's considered 
previously, since there is a privileged x direction. Another simplification would be to 
consider a coupling of the large spin j with a large collection of independent systems, as 
done in [Strunz et al., 2002 . However, here again, to do so while preserving the SU(2) 



invariance of the ensemble of random interaction Hamiltonian s is not that easy and does 
not lead to models simple enough to be discussed here. 

We have not studied the case where the internal dynamics for the spin is non-trivial 
[Hg 7^ 0). This has been done analytically for the spin j = 1/2 case in |Lebowitz and Pastur, 2004 



Lebowitz et al., 2007| , and they have obtained explicit expressions (of course more com- 



plex than our results for j = 1/2). It should not be difficult to extend their results for 
our class of interaction Hamiltonian in the specific j = 1/2 case (then only the / = 1 
sector contributes). However the problem becomes increasingly complicated as the spin 
j increases, and no general solution valid for generic spin j has been obtained yet. This 
is clearly an interesting problem, relevant for studying the interplay between decoherence 
and dissipation. If the internal spin dynamics is slow, i.e. when H$ is small, it should be 
possible to study the dynamics by standard approximation methods (short time expan- 
sion, TCL). 

We have studied the matrix model in the planar large iV limit, where iV corre- 
sponds to the size (number of states) of the external system. It should be interest- 
ing to study the case of finite but large N, as has been done already for j = 1/2 in 
|Esposito and Gaspard, 2003 . Especially interesting should be the case where N (the 



dimension of the external system) is of the same order than j (the spin). 

In our model we have considered a random M x M matrix model of the GUE type, but 
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such that the distribution probability is not invariant under the whole U(M) group, but 
only under the action of a subgroup G acting on the space (in our case G =SU(2) xU(iV) 
with M = (2j + 1)N). This class on models belongs of course to the very general class 
of block random matrix models, but in our case the presence of the symmetry group G 
and of its spin representation adds a lot of structure and of simplifications. It should be 
interesting (at least mathematically) to study these kind of models for more general rep- 
resentations of the group G =SU(2) xU(iV), and for general groups and representations. 

Finally, the model presented here is very idealized and is mostly of pedagogical and 
mathematical value. But it should be a first step in studying general but more realistic 
models of spin decoherence, in particular for systems where the classical macroscopic 
degrees of freedom do not come from a per se large quantum spin, but emerges from the 
interaction between many small quantum spins. 
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The function Z(l) for various coupling distributions 
A(Z) and spin j 




Figure 16: Z(l) for various Iq and j, all A(Z)'s are equal 
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coupling distribution A(l) = {0, 1} 
total spin j = (1, 2,4, 8, 16, 32, 64, 128| from blue to red 



coupling distribution A(l) ={0,1,1} 
total spin j = {1, 2, 4, 8, 16, 32, 64, 128, 256| from blue to red 




coupling distribution A(l) = {0, 1, 1, 1} 
total spin j = {3, 6, 12, 24, 48, 96, 192, 384} from blue to red 



coupling distribution A(l) = {0, 1, 1, 1, 1} 
total spin j = {2, 4, 8, 16, 32, 64, 128, 256, 512} from blue to red 




coupling distribution A(l) = {0, 1, 1, 1, 1, 1} 
total spin j = {5, 10, 20, 40, 80, 160, 320, 640| from blue to red 



coupling distribution A(l) = {0, 1, 1, 1, 1, 1, 1} 
total spin j = {3, 6, 12, 24, 48, 96, 192} from blue to red 




its 1.0 



Figure 17: Z(l) for various Iq and j, all A(7)'s are equal but A(0) = 
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coupling distribution A(l) = (0, 1} coupling distribution A(l) = {0, 0, 0, 1} 

total spin j = (1, 2,4, 8, 16, 32, 64, 1281 from blue to red total spin j = {3, 6, 12, 24, 48, 96, 192, 384| from blue to red 




-1.0 - ' -1.0 



Figure 18: Z(l) for various /q and j, one odd I only 
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Figure 19: Z(l) for various Iq and j, only odd Z's , all equal 
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coupling distribution A(l) = (0, 0, 0, 1, 1} 
total spin j ={2,4, 8, 16,32,64, 128,256,5121 from blue to red 



coupling distribution A(l) = (4, 0, 0, 1, 1} 
total spin j = {2, 4, 8, 16, 32, 64, 128, 256, 512) from blue to red 




coupling distribution A{1) = {16, 0, 0, 1, 1} 
total spin j = {2, 4, 8, 16, 32, 64, 128, 256, 512) from blue to red 



coupling distribution A{1) = {64, 0, 0, I, 1} 
total spin j = |2, 4, 8, 16, 32, 64, 128, 256, 512) from blue to red 




0.2 0.4 0.6 0.8 1.0 



02 0.4 0.6 0.8 1.0 



Figure 21: A(0) becomes large, the others are fixed 



coupling distribution A(l) = {0.841765, 0.291879,0.670747,0.153111,0.190376) coupling distribution A(l) = {0.367858,0.61796, 0.654416, 0.0657948,0.916895) 
total spin j = {2, 4, 8, 16, 32, 64, 128, 256, 512) from blue to red total spin j = |2, 4. 8, 16, 32, 64, 128, 256, 512) from blue to red 

zi/i m 




Figure 22: Random A(7)'s 
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